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Abstract 


This report is a self-contained and comprehensive review of the physics of 
propagating pulses of high-intensity high-energy particle beams in pre-existing 
or self-generated plasmas. The first part is a systematic overview of the theory 
pertaining to propagation, plasma self-generation, energy/current-losses, and sta- 
bility of such pulses. The second part reviews the major full-scale propagation 
experiments which have been carried out, in atmospheric and outer-space plas- 
mas, to assess the validity of theoretical models. It is found that the data available 
on these experiments demonstrate that range and stability are in agreement with 
theory. In particular, stable self-pinched propagation of high-current beams in the 
atmosphere is possible over distances on the order of several Nordsieck lengths. 
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Chapter 1 


Introduction 


This report deals in a comprehensive manner with theoretical and experimental 
plasma-physics and accelerator-physics research which have been actively fol- 
lowed or done by the author over the past twenty-five years, that is starting from 
approximately the time of the beginning of the construction of the “Advanced Test 
Accelerator” (ATA) at the Lawrence Livermore National Laboratory, and of the 
creation of the “Accelerator Technology” division (AT) at the Los Alamos National 
Laboratory. 


In order to understand the relevance of the key experiments, which started to 
give significant data in the mid-1980s, it is important to understand the theory 
underlying the numerous plasma-physics effects at work when a high-intensity 
beam-pulse of particles propagates in a background gas or through the atmosphere. 
Since there is no published text book or monograph covering this subject in a 
systematic manner, Chapters 2 to 7 attempt to make a synthesis of numerous 
published articles and many informal reports which deal with one or another 
aspect of this theory. As shown by the bibliography, this meant studying many 
papers published over the past fifty years, often dealing only indirectly with the 
subject, in order to extract the pertinent information necessary to make a consistent 
theoretical model. 


In the second part of the report, Chapters 8 to 10, the focus is on the acceler- 
ator facilities and the beam propagation experiments which are significant for the 
purpose of establishing the feasibility of generating suitable high-intensity high- 
energy particle beams, and of propagating them in outer-space plasmas or through 
the atmosphere. The difficulty, here, is that the openly available data is more of 
a qualitative than quantitative nature, which is precisely why a thorough under- 
standing of the plasma-physics pertinent to these experiments is so important. In 


these chapters, which deal with technologies at the frontier of the state-of-the-art, 
an effort is made to refer to the implications of the most advanced theoretical 
ideas and technologies, in order to show how much the possible future engineering 
development of high-power particle beam generation/propagation technology still 
depends on ongoing research. 


While we said that there appears to be no published text book or monograph 
covering the subject of this report in a systematic manner,’ there is a growing 
number of excellent books available on the physics of charged particle beams and 
their applications. The books by J.D. Lawson [I], R.B. Miller [2], S. Humphries 
Jr. BIA, and M. Reiser [5], are possibly the most useful in the context of the 
present report. 


As the information and research summarized in this report extend over so many 
years, there are many people to thank for their direct and indirect contributions 
to it. While I cannot mention all of them, I wish in particular to thank my 
former colleagues at CERN (where this work started): Claude Bovey, Steve Geer, 
Peter Jenni, Pierre Lefèvre, Claude Metzger, Dieter Moehl, Emilio Picasso, Peter 
Sonderegger, Daniel Treille, and Horst Wachsmuth; as well as Frank Barnaby and 
Bhupendra Jasani at SIPRI (where part of this report was written); and last but not 
least, Jean-Pierre Hurni at ISRI. 


'One exception appears to be the lecture notes prepared by Prof. K.E. Woehler, Department of 
Physics, Naval Postgraduate School, Monterey, CA, for his course PH 4959 — Physics of directed 
energy weapons: Part I, Particle beam weapons (March 1981) 120 pp; Part II, Particle accelerators 
(March 1981) 60 pp. However, the level and scope of these lectures are more elementary and less 
comprehensive than those of the present report. 


Chapter 2 


Some preliminary definitions 


A particle beam pulse may be thought of as an ensemble of moving particles whose 
trajectories constitute a “bundle.” The diameter of this bundle is small compared 
with its length, and the trajectories generally make a small angle with the “axis.” 
The complete description of the evolution of such a system of interacting particles, 
especially if they propagate through a gas or plasma, is in general very complicated. 
However, in many cases, the beam pulse can be characterized statistically by the 
RMS (i.e., “root mean squared”) values of its radius, length, angular spread, 
energy spread, etc. A good description is then provided by the so-called envelope 
equations giving the RMS radius (or length) of the pulse as a function of time or 
propagation distance. 


A most important concept specifically related to particle beams is that of 
emittance, which can be considered as a measure of the disorder in the motion of 
the particles relative to the average motion of the beam [6]. In the general case 
it is necessary to distinguish between the transverse and longitudinal emittances, 
i.e., 2 and ,, which will be defined in the first section of the next chapter. Since 
an emittance is basically the product of the extent of a spatial distribution along a 
given direction by an angular spread, i.e., a length times an angle, it is generally 
measured in units of “m rad.” 


A concept that is related to beam emittance is that of beam brightness, i.e., 
current divided by angular beam spread, which means that low-emittance implies 
high-brightness, and vice versa. Both concepts are often used interchangeably to 
qualify a high-power directed beam, but we will use only the concept of emittance 
in the theoretical sections of this report. 


In order to simplify the calculations, the usual treatment generally assumes 
that v. and %, the RMS values of the random components of the transverse (or 


perpendicular) and longitudinal (or parallel) velocities, are small compared with 
the mean longitudinal drift velocity v= c: 


2o d; 22: (2.1) 


This is the paraxial approximation and the particle trajectories deviate only slightly 
from parallel straight lines. In such a model, the beam particle momentum p, 
total energy W , and kinetic energy K are slowly changing parameters with the 
longitudinal distance. With = 1= 1 2 the Lorentz factor, 


p= mc; W= mé;K=( 1)mé: (2.2) 


In the case of charged particle beams, we will see that the paraxial approxima- 
tion is equivalent to the statement that the effective beam current I, generating the 
electromagnetic self-fields is small compared with the Alfvén current I, defined 
as [7] : 
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where q= Z Fjis the electric charge of the beam particles. 


In the discussion of problems like beam-plasma intęgaction and stability, the 
most convenient radial scale is not the RMS radius a= r°, but the scale radius 
a defined such that 


Js (0) = — (2.4) 
where Jg (r) is the current density and I, the total beam current. 


We will use a cylindrical coordinate system with radial distance r = j x? + y?, 
azimuthal angle and longitudinal distance z. In the beam frame, we will replace 
the longitudinal variable z by a time coordinate which is zero at the beam head 
and equalto __, the pulse duration, at the tail-end. 


In the following sections we will first review the envelope equations for neutral 
and charged particle beams. We will assume homogeneous and constant back- 
ground conditions, and neglect transient effects and instabilities. Second, we will 
examine the transient effects at the head of a charged beam when it is fired into 
a initially neutral gas such as the atmosphere which is turned into a plasma by 
the beam. Finally, we will examine the main possible instabilities affecting the 
propagation of such a beam. 


Chapter 3 


Particle beam propagation in a 
vacuum or a tenuous plasma 


3.1 Neutral beam in a vacuum 


In the absence of collisions between the particles, the individual trajectories of 
neutral particles propagating in a vacuum are all straight lines. The envelope 
equations can then directly be derived from kinematics [8]. For a beam pulse with 
axial symmetry, the RMS radius a, and the RMS half-length * are given by: 


2 
k 
ae: (3.1) 


In these envelope equations, the primes denote derivation with respect to 
the longitudinal coordinate z. The constant p and , are the transverse and 
longitudinal emittances which characterize the random distribution of the particles 
in the beam. Specifically 


2 a’ WA )? )? P 2 2 pea (%)? cy g 


2 rok 


(3.2) 


v? y 
As 2 and , are constant of motion, they can be measured most conveniently at 
the point where the beam envelope forms a waist. At such a point a°= Oor “= 0, 
and 

Y W 

2 W ° 
where, in the paraxial limit, ~is the RMS angular spread and Ww =w the RMS 
energy spread of the beam. The emittances > and , are thus characteristic of the 


2 = a~ 7 k7 (3.3) 
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beam quality: a small transverse emittance corresponds to a well collimated small 
angular divergence beam, and a small longitudinal emittance to a short pulse with 
a small energy spread. 


For possible applications for space-based beam weapon systems, the solutions 
of B.I) have one essential feature: for a given emittance, and at a fixed distance z 
from the accelerator generating the beam, the minimum beam spot size is inversely 
proportional to its initial radius ao: 


Qn in = Z— (3.4) 
&o 

This minimum is achieved for a} = a =z and shows that the critical parameter 
in focusing the beam is its emittance and not just its angular spread or radius. 
Similarly, the longitudinal emittance determines the minimum duration of beam 
pulse, and is thus an essential parameter for compressing the beam into short pulses 
at the target. According to 8.4), in order to focus a neutral beam into 1 m radius 
at 1000 km, one would, for example, need a beam with an initial radius of 20 cm 
and an emittance of 2 10’ m rad. An accelerator with such an emittance, and 


a focusing system for such a beam, are currently under development at the Los 
Alamos Scientific Laboratory in the USA [9]. 


3.2 Charged beam in a vacuum 


When a charged particle beam is launched into a vacuum, the beam will tend 
to spread apart because of its emittance and the effect of the Coulomb repulsion 
between like charges, which is particularly strong if the beam energy is low. 
However, for relativistic beams, this is not the main problem with propagation 
in a complete vacuum: when an unneutralized beam leaves the accelerator, the 
accelerator becomes charged with the opposite sign, pulling the beam particles 
back and causing them to decelerate. Let us examine the importance of these 
effects. 


In the paraxial limit the radial electrostatic force on a particle within a beam 
at a distance r from the axis is 


W i (£) 
r 


7 (3.5) 


where I, (r) is the total current flowing within the radius r. The outward-directed 
force is partially compensated for by the inward-directed pinch force due to the 
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azimuthal magnetic field generated by the beam current this radius: 


1 2W ke) 
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(3.6) 


As ? < 1, the net force Fẹ + F, is outward-directed, and causes the beam to 
spread apart. 


The radial envelope equation combining the effect of the emittance and those 
of the electro-magnetic self-field is [0]: 


2 
ns Seg (3.7) 


In order to estimate the effect of the electromagnetic spreading term, one can first 
neglect the emittance in (3.7). The distance over which a beam has to propagate 
in order for it radius to double is the approximately given by [LL]: 


pP 
Z2 & Is =Ia (3.8) 


With an initial radius of a= 20cm, a kinetic energy of 500 MeV and a current of 
100 mA, this distance would be about 5000 km for electrons, but only 5 km for a 
proton beam with the same characteristics. But this is estimated from radial beam 
divergence alone. In fact, as there is no return current the particle will also be 
slowed down under the action of the longitudinal restoring force due to charging 
up of the spacecraft launching the beam. This effect puts a limit to the range of 
charged particles in a vacuum. At relativistic energies [12], 


r 


Tp 
2 ae 3.9 
Žana & 2, (3.9) 


and both proton and electron beams of the above example would have a range of 
less than 1 km. This calculation alone would tend to rule out the use of charged 
particle streams as possible beam weapons in a vacuum. However, as will be 


examined now, if such beams were injected into outer space, which is in fact a 
dilute plasma and not a vacuum, the situation can be very different. 
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Chapter 4 


Particle beam propagation in a 
dense gas or plasma 


4.1 Charged beam in a plasma : the Bennett pinch 


Charged particle beams for use as directed energy weapons are in general injected 
either into the atmosphere for ground- and aircraft-based systems, or into a plasma 
for space-based systems. In the case of atmospheric systems, the beam will enter an 
initially neutral atmosphere and, by ionizing the air, turn it into a plasma along the 
beam path. In the case of outer-space systems, the plasma will be the ionosphere 
for near earth orbiting systems, or the interstellar environment. In both cases, the 
transient phenomena occurring at the head of a beam pulse are complicated. We 
will thus concentrate first on infinitely long beams in the paraxial approximation. 
Similarly, we will assume that the plasma can be sufficiently well described for 
our purpose by a scalar conductivity given by 


s2 (4.1) 


where is the collision frequency, m e the rest mass and n, the number density of 
the conductivity electrons. 


In order to appreciate the relative importance of the plasma background for 
beam propagation, it is sufficient to compare the plasma density ne to the beam 
particle density n,. On the beam axis, according to B1), 


1 
Np (0) = as (4.2) 


For example, in the case of a relativistic beam with I, = 100mA and a= 20 
cm,n,= 1:6 10cm °. In comparison, in the ionosphere between an altitude 
of 100 to 2000 km, the electron number density ne is on the same order or larger. 
Therefore, plasma effects cannot be ignored. 


The first effect of a plasma background is that, because of quasi-neutrality, 
the excess charge locally introduced by the passage of the beam will tend to be 
neutralized on a time scale į set by the plasma conductivity: 


so (4.3) 


As a result, the space charged force will be reduced to a value 


Fe= il £)— (4.4) 


where feis the space charge neutralization fraction. 


The actual physical processes involved at the microscopic level are different 
depending upon the sign of the electric charge of the beam particles. For an 
electron beam, the electric field will quickly expel the plasma electrons and the 
neutralization will be provided by the positive ions left within the beam. Assuming 
the ions not to be able to move significantly during the passage of the beam pulse, 
the maximum charge neutralization fraction will be fe = 7py=n,jin cases where 
ni < ny. In the case of a proton beam, electrons from the plasma surrounding 
the beam will be attracted into the beam region. The neutralizing fraction can be 
1 even when ne < n,. The fact that the charge neutralizing particles have very 
different masses means that the transient phenomena and stability conditions can 
be different in the two beams. However, in dense plasma, i.e., when ne n> np, 
most properties will be similar. In particular, fe will essentially be one for beam 
pulses of duration longer than «and radius smaller than œ e. 


The second main effect of a plasma background is that return current can be 
established. This current is driven by E z, the longitudinal electric field induced 
by the beam. From Faraday’s law of induction 


1 @ 


E, = — 
oZ At 


CI) (4.5) 


where I; is the effective current driving the electromagnetic self-fields and L a 
dimensionless inductance. For axially symmetric beams 


1 loa 
L grlt a (4.6) 
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where bis the maximum radius out to which the plasma background is significantly 
affected by the beam. For instance, in beam-generated plasmas, b is normally 
determined by the extent of induced breakdown around the beam head. Typically, 
ba  10,1e.,L 1. At the boundary of this region the conductivity becomes too 
small to insure quasineutrality. The charge imbalance from the beam is conducted 
to this surface, which is thus the path along which the beam current not neutralized 
by the plasma current is returned to the accelerator. For Ohm’s law, J = E, the 
longitudinal component of the plasma current is then 


@ 
b= a E: = nate) (4.7) 
where 7 
a 
m = —— 4.8 
g (4.8) 


is by definition the magnetic diffusion time of the plasma. From it turns out 
that the plasma current will in effect decrease the beam current and thus reduce 
the magnetic pinch force. This is usually written as 


(4.9) 


where fn = jp =Ñ jis the magnetic neutralization fraction. Because E, is a 
function of the beam current and shape variations, the plasma current, and thus 
fn , will be largest at the beam head. Within the pulse, fw will decrease to zero on 
a time scale set by the magnetic diffusion time. 


Assuming feand fm to be independent of r, the envelope equation including 
the plasma effect is then [10] 


(4.10) 


where I, is the effective current associated with the net electromagnetic force on 
the beam: 


1 
E=2 a f) =a f) : (4.11) 
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The effective current contains charge imbalance as well as true beam and plasma 
currents. It can have both signs, and the forces can either tend to separate or, on 
the contrary, to pinch the beam. When J, is positive, a stationary solution with 
a= a= Qis possible. This is the Bennett pinch solution, and in that case 
gives [I3]: r 


a=? > (4.12) 


Figure 4.1: The pinch effect. The photograph, taken in the early 1950s at the 
Argonne cyclotron (near Chicago), shows the glow produced when sending the 
full deuteron beam into the atmosphere. Because of ionization the air near the 
beam is turned into a plasma which keeps the beam from expanding radially 
under the effect of Coulomb repulsion between like-charged particles: This is the 
“Bennett pinch effect,” first described by Willard H. Bennet in 1934 to explain 
focusing effects and breakdown in the residual gas of high-voltage electronic 
tubes, and later applied to the propagation of interstellar and interplanetary self- 
focussed beams of particles, such as proton streams traveling from the sun towards 
the earth. As the beam loses energy and intensity because of interactions with 
atmospheric nitrogen and oxygen nuclei, plasma generation becomes less efficient 
and the beam progressively expands: This is the “Nordsieck effect,’ after the 
name of Arnold Nordsieck who is generally credited for having first explained 
this expansion. Ultimately, when the plasma effects become to weak to pinch the 
beam, it breaks-up. This happens at a propagation distance on the order of the 
so-called “Nordsieck length.” 


From (4.11) we see that the Bennett pinch radius (4.12) is minimum for f = 1 
and fw = 0. In that case I, = I, the pinch force is maximum and the beam 
is fully pinched. The Bennett pinch solution exists, however, only in the paraxial 
limit. This can be seen from and which imply 


aL (4.13) 


In a Bennett pinch, the beam particles perform harmonic motion around the 
beam axis. The angular frequency of the rotation is a function of r and is called 
the betatron frequency ! . By equating the net force Fe + Fm to the centrifugal 
force m !?rand by averaging over the beam current density one finds from 


15 


and that the mean azimuthal velocity is 


I 
yal ya (4.14) 
Ta 
which is independent of the beam profile. In the general case, the betatron 
frequency ! * is distributed between zero and a maximum 


17 = !70)= 2 75 = (4.15) 
which is also independent of the beam profile. In the special case of a beam with 
a constant current density out to a radius a, the betatron frequency is constant and 


equal to (4.15). 


For propagating self-pinched beams, the most natural equilibrium density 
profiles are those corresponding to a Maxwellian (i.e., Gaussian) transverse energy 
distribution [14]. 


Possible equilibria include filamentary current flows along the axis, but the 
simplest realistic example is the so-called Bennett distribution [I3] 


Jp (Ct) = — 1+ — : (4.16) 


The RMS radius of this distribution diverges logarithmically. However, both 
theory and experiment that indicate a Bennett profile are not valid for 
r a. In practice, the current profile is often considered to be truncated at 
r= 2a. This yields a = 1:006 a. 


In a beam with the Bennett profile (4.16), the betatron frequency !?° is is 
uniformly distributed between zero and the maximum given by (4.15). 


The significance of the Bennett pinch is that when a beam is injected in a gas 
or plasma, sufficiently dense to suppress the effect of the space charge repulsion, 
the beam may pinch down to a minimum constant radius and propagate over large 
distances. In particular, because the plasma background provides a means for 
carrying the return current, the range is no more strictly limited as ii the case of 
charged beams in vacuum. The limit to the propagation distance will now be set 
by scattering, energy loss, instabilities, etc. 


The fact that the net charge transported by a beam pulse traveling through a 
plasma is equal to the charge of the beam itself, even when the beam is launched 
from a ground plane, is not an obvious result. It is, however, correct, even when 
the plasma is generated within the pulser by beam-gas interaction [I7]. 
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The Bennett pinch is important in many area of science and technology. In 
particular, it is important for understanding interplanetary particle streams 
and studying the ionosphere with beams launched from rockets or the space 
shuttle [T9]. It has many applications in thermonuclear fusion research and has 
been envisaged as a means for accumulating high energy electrons in large rings 
in outer space. Such rings could be used to store energy and subsequently gen- 
erate synchrotron radiation or free electron laser optical beams [20]. Finally, for 
endo-atmospheric beam weapon, the pinch effect provides the means for radially 
confining charged particle beams. 


The Bennett pinch condition I, > O requires fẹ > 1= * when fa = 0. For 
high energy beams, i.e., = 10 16, this condition is easily satisfied, even in very 
low density plasmas. This allows the transport of beams with current densities 
higher than the plasma density by a factor ?. 


The practicality of this method for stable propagation of beams over large 
distances in a low pressure atmosphere was demonstrated for electrons in 1972 
BI]. For beams such as those considered at the end of the previous paragraph, 
the pinch effect would provide an ultimate theoretical range of many thousands of 
kilometers. However, the practical use of such beams in outer space will, in fact, 
be limited by the effect on them of Earth’s magnetic field. On the other hand, in 
a full density atmosphere, the possibility of using the pinch effect to propagate 
a charged-particle beam has still to be experimentally demonstrated. This is the 
goal of various large scale experiments currently underway in several laboratories 


(22) 23}. 


4.2 Effect of magnetic fields : neutral plasma beams 


The region above the atmosphere in which an orbiting charged particle beam 
weapon might be deployed contains plasma and magnetic fields of both solar and 
terrestrial origin. The plasma may enable electron or proton beams to pinch and 
propagate over large distances, but the magnetic field would strongly affect the 
beam trajectory in most cases. 


In a homogeneous magnetic field B the radius of curvature of a particle of 
charge qis 
P 
= — (4.17) 
QB 
hence, for an electron or proton with momentum p = 1 GeV/c and in a field of 
0.3 Gauss typical of Earth’s geomagnetic field R would be 110 km. In the 


absence of any other effect, such a deflection would be a small correction for 
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an endo-atmospheric system with a range of a few kilometers. However, for an 
outer-space system with a required range of several thousands of kilometers, the 
beam would spiral much as it did in the low energy beam experiment performed 
on the space shuttle [19]. 


For beam without current neutralization (i.e., fn = 0), the effect of an external 
magnetic field is essentially the same as that on the beam particles taken individ- 
ually. In that case, applies directly to the beam as a whole, which has been 
verified experimentally [24]. On the other hand, if fa 6 0, the overall effect of 
an external magnetic field on a beam would be smaller, and would vanish in the 
extreme case of a beam with fully neutralized current. 


This phenomenon is most easily understood in the ideal case of a neutral 
plasma beam consisting of co-moving electrons and ions: electrons and ions are 
deflected in opposite directions and create a polarization electric field. This field 
exactly cancels the force of the external magnetic field and the beam continues 
undeflected, loosing, however, particles from the polarization layers [25]. This 
effect has been recently investigated experimentally for a charge-neutral beam 
consisting of electrons and protons [26]. For such a beam, the effective current 
is obviously zero. The envelope equation reduces then to that of 
neutral beam and, in the absence of pinch effect, the neutral plasma beam would 
expand through the effect of its emittance [27]. 


A neutral plasma beam, also called plasmoid beam, would thus be an alternative 
to a neutral particle beam for outer space systems. The main problem is to generate 
such a plasmoid beam with a sufficiently small emittance for the range to be useful, 
and to overcome the problem of the beam losses at the boundary layers. 


In an atmosphere with a pressure of about 1 Torr, stable propagation conditions 
with fe = fa = 1 were discovered in early electron beam experiments, first in the 
United States [28]29], and later in the Soviet Union [12]. In these experiments, the 
current neutralization was provided by slowly counterstreaming plasma electrons. 
The fact that such a configuration is stable and unaffected by an external magnetic 
field would make it attractive, in principle, for propagating a beam in a reduced 
density channel. However, the beam is not pinched, and would thus quickly spread 
apart because of collisions with the gas molecules. 
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4.3 Charged beam in a dense gas or plasma : Nord- 
sieck equation 


When a beam propagates in a dense gas or plasma, the collisions with the electrons 
and nuclei result in a loss of energy by the beam current when particles can 
effectively be removed from the beam by interactions with nuclei. Furthermore, 
by increasing the angular and energy spread of the beam, the collisions lead to a 
continuous increase in both transverse and longitudinal emittance. In the paraxial 
approximation the full radial envelope equation taking energy loss and scattering 
into account is as follows [O {15}: 


(4.18) 


This result is obtained from kinetic theory which shows that a self-pinched 
beam subject to gas scattering evolves to a state in which its current density takes 
the form of the Bennett distribution (4.16), which is a similarity solution of the 
Fokker-Planck equation witha aforr O ). 


In all parameters are, in general, a function of the propagation distance z. 
The right hand side corresponds to the increase in transverse emittance because of 
Coulomb multiple scattering in the background gas. In agreement with Liouville’s 
theorem [fi], the emittance enters this formula through the product pe which is 
the relativistic invariant conserved emittance for a beam of varying energy. The 
multiple scattering angle is given by Rossi’s formula p.67] 


Es 2 AZ 
cp Xo 


d= (4.19) 


pP F er 
where E= 4= mœ 212 MeV and X) = 4 x Na =A is the radiation 


e 


length of the medium, i.e., X 9 300:5 m for air at standard temperature and 
pressure (STP).'‘? 


The total energy of the beam particles as a function of zis given by the equation 


dw W 
— = — +c, : 4.2 
= SO 7 + (4.20) 


The first term is the Bethe’s stopping power and corresponds to energy losses 
by ionization and excitation of the gas molecules [B0]. This loss is a slowly 
! = 1=137 is the fine-structure constant, and rẹ = e?=m eœ the classical electron radius. 

Rossi’s formula provides a convenient first approximation for analytical calculations of the 
kind done in this section. For more precise calculations, especially for beams of relatively low- 
momentum particles, it is better to use Moliére’s theory as formulated by Bethe [BT]. 
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varying function of energy and S equal 2 to 3 keV/cm in air at STP for 1 GeV 
single-charged particles. The second term is present only for electron beams and 
corresponds to bremsstrahlung radiation losses [30]. (Bremsstrahlung emission is 
negligible for particles of mass heavier than electrons). The third term corresponds 
to the ohmic losses in the induced longitudinal electric field. These losses which 
are associated with the charge and current neutralization process are concentrated 
in the head of the beam where they contribute to heating the plasma electrons. 


The beam current decreases because of collisions with nuclei. In the case of 
proton beams, both inelastic and elastic nuclear collisions will effectively remove 
particles from the beam. Therefore, in that case, the effective beam current 
will vary with z as 

Ir @) = Tk 0) exp( z=Xn) (4.21) 


where X n is the nuclear collision length, about 500 m in air at STP for high 
energy proton beams. In the case of electron beams, the pronounced statistical 
character of bremsstrahlung radiation losses will result in a wide energy spread for 
propagations over distances of the order of one radiation length. The implication 
of such a large increase in longitudinal emittance is that particles with energies 
less than half the mean beam energy will in fact “evaporate” from the beam, thus 
leaving behind a reduced current beam [32]. A calculation based on Bethe and 
Heitler’s theory of straggling in bremsstrahlung emission [30] indicates that this 
effect would account for a 30% loss in beam current for a high energy electron 
beam propagating over one radiation length. 


The solution of has generally to be found numerically. However, con- 
siderable insight can be gained in the quasistatic limit which is taken from the 
Bennett pinch condition where a a. In this limit, $ = a°= 0, so that 
becomes the so-called Nordsieck equation (10) (15): 


Py 1 
oie. a)y= Æ (4.22) 


where Pg = Jy (z)p(z)ee and Py = H2=e = (rn) 4 = mc 15 
10’ W = 15 TW is a constant: the Nordsieck power. The remarkable property 
of this approximate envelope equation is that the radial expansion of a relativistic 
beam in a gas is a function of only two variables: the radiation length X ọ which 
characterizes the medium, and the effective power Pg, which characterizes the 
beam. 


The Nordsieck equation can be solved explicitly for several cases of inter- 
est. For instance, for high energy proton beams, the dominant beam power loss 
comes from the decrease in beam current due to nuclear interactions: P (z) 


20 


kosik CE Yc a BY 


400 4000 40000 (m) 


Figure 4.2: Beam expansion for various particles. The radial expansion of a 
beam of intensity IZ = 10 kA and momentum Pg = 10 GeV/c, calculated by 
integrating numerically the full radial envelope equation taking all energy losses 
and scattering effects into account, is plotted as a function of propagation distance 
in air at STP. The efolding range, also called the Nordsieck length, defined as 
the distance at which the beam radius has expanded by a factor of e 2:718, is 
of about about 400, 1000, or 4000 meter for an electron, proton, or muon beam, 
respectively. 


21 


(m) 
4 500 


4000 


500 


0.1 3j j0 50 


Figure 4.3: Beam range versus momentum for various particles. The effective 
range, defined as the distance at which the beam radius has expanded by a factor of 
2e 5,is plotted as a function a beam momentum in GeV/c. For a beam intensity 
Is = 10 kA the range in air at STP of a single 1 GeV/c electron, proton, or muon 
beam pulse is of about 200, 500, or 800 m, respectively. For a 10 GeV/c pulse the 
effective range increases to about 500 or 1200 m for an electron or proton beam, 
and about 6000 m for a muon beam (see Figure 2). 
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P,yexp( z=X,). Then 


1 Z y Enžn 


a= a a Á= 
SeS Aa Poki 


exp (Z=X,) 1 (4.23) 


On the other hand, for high energy electron beams, the dominant effect is energy 
loss by bremsstrahlungs: Pz (z) P exp( z=Xo). Thus 


a= a exp- — + — exp(z=Xo) 1 (4.24) 


In both cases, the Nordsieck power Py plays an important role. For Py < Py, the 
radial expansion is very fast and the beam cannot propagate over sizable distances. 


For Po > Py , however, the efolding range, also called the Nordsieck length, 
defined as the distance over which the beam radius expands by a factor of e 
2:718, becomes independent of initial beam power and is about 2X„ 1000m for 
protons and 2X 5 600 m for electrons in air at STP. The implication is that for a 
beam to propagate, its total effective power has to be of the order of Py = 15 TW, 
but any further increase in beam power would not substantially increase its range. 


In the low-energy limit, i.e., Py < Py, the beam will propagate less than one 
radiation or nuclear interaction length. One can therefore assume that z X for 
electrons,andz X „ for protons, to find that in first approximation the Nordsieck 
length for both an electron or a proton beam is given by 


Po 
———_— g : 
Zy T Fi (4.25) 
For example, the maximum power of the 10 kA, 50 MeV, ATA beam is Po = 
0:5 TW. Since the Nordsieck power is Py = 15 TW, and the radiation length 
Xo 300 m in air at STP, the typical propagation length of a single ATA beam 
pulse should be on the order of zy 20m. 


In order to check the validity of these conclusions, we have solved for 
10 GeV, 10 kA beams of electrons, protons and muons in air at STP, modeling 
all physical effects such as multiple scattering and energy losses as precisely as 
possible. (See Figure 2.) This corresponds to an initial beam power Po = 6:6Py. 
The range of the proton beam is about 1 km, as expected. But the range of the 
electron beam is only 400 m. This discrepancy is the result of the approximation 
made in deriving the Nordsieck equation and of neglecting straggling. The 
muon beam, which would be much more difficult to produce in practice than an 
electron or proton beam, is included for comparison. As bremsstrahlungs and 
nuclear interactions are negligible for these particles, the muon beam range is of 
course much larger. (See Figures 2 and 3.) The range of proton beams in air is 
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thus strictly limited by the nuclear interaction length, and electron beams by the 
radiation length. However, as these lengths depend directly upon the atmospheric 
density, the corresponding ranges will increase in proportion to the decrease in 
atmospheric density. 


Finally, in order to compare the beam expansion theory given in this section to 
the data, and therefore to validate the extrapolations given in Figures 2 and 3, we 
have simulated and found good agreement with the beam expansion measurements 
given in reference [16]. In these measurements, the ASTRON induction linear 
accelerator of the Lawrence Livermore National Laboratory was used to propagate 
a 250 ns, 0.85 kA, 5 MeV electron beam pulse over a distance up to 16 meters in 
reduced density nitrogen. This agreement can be also be used to predict with some 
confidence the propagation characteristics of the 10 kA, 50 MeV, ATA beam. It is 
then found that one e- folding beam-radius increase corresponds to a propagation 
distance of about 20 m in a sea-level air, in excellent agreement with the analytical 
estimate given by equation (4.25), and that the beam radius increases to about 10 
times its initial value after propagating a distance of about 35 m. 
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Chapter 5 


Injection of a high-power beam into 
the atmosphere 


5.1 Plasma generation by a particle beam 


In a endo-atmospheric system, except when the beam is launched into the iono- 
sphere or a pre-formed plasma channel, the beam will, in general, be injected into 
initially un-ionized air. At the very head of a beam pulse, there is thus no back- 
ground plasma and therefore, both fẹ and fa are zero. The beam head, therefore, 
expands at a rate governed by the net radial electromagnetic force, emittance and 
scattering. However, as the beam ionizes the air, a plasma of increasing density 
builds up and, as plasma currents start flowing, the space charge gets progressively 
neutralized. When *(f, *  £)= 1, Ik = O, the net electromagnetic force 
changes sign, and the beam stops expanding and starts pinching down. During the 
self-pinching of the beam, both the radius a and the effective beam current I, are 
rapidly varying, and this happens while the beam current Ip, which in practice 
has a finite rise-time, increases. As a result, a strong electric field is induced. 
This field, given by (4.5), accelerates the plasma electrons, which, as they gain 
sufficient energy, start ionizing the gas as well. The rate of neutralization of the 
beam thus increases, until it levels off as the decreasing beam radius approaches an 
equilibrium radius given by (4.12) with fe 1. At this point, the induced electric 
field becomes very small and, the plasma current I, , after having gone through 
a maximum, starts decreasing according to (4.7). Finally, for a sufficiently long 
pulse, after the plasma current has completely decayed, the beam becomes fully 
pinched. 


As a result of this process, the beam pulse takes on a characteristic “trumpet” 
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Figure 5.1: Propagation of a pinched beam. The diagram shows a propagating 
beam indicating the pinching, necking, and freely expanding regions. In the high- 
conductivity pinched (also called “Nordsiecking’”’) region the particles perform 
betatron orbits with a slowly increasing radius. In the neck region the beam 
expands and undergoes erosion as particles are lost as a result of weakening pinch 
forces. Finally, in the low-conductivity free expansion region where the plasma 
effects are not sufficient to pinch the beam, the particles move away from each 
other because of Coulomb repulsion. 
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Figure 5.2: Coordinate system used in describing the propagating beam. The 
variable z is the distance from the accelerator to a point within a beam pulse, while 
the variable is the time such that the distance from this point to the beam head 
is equal to v, , where w is the longitudinal velocity of the beam’s particles. 


shape, and can be divided into four distinct regions; the expanding beam head, the 
neck region in which E, and f, are maximum, the body in which fe 1 and 
Nordsieck’s equation becomes a good approximation of (4.18), and finally the tail 
where fe = land fm 0. (See Figure 4.) 


The description of this process requires an equation for the plasma electron 
density ne and temperature Te so that the conductivity s, and thus .and m can 
be calculated. Furthermore, in order that all quantities can be expressed self- 
consistently an equation relating the various fields and current is needed. This last 
equation is called the circuit equation. These equations, which describe quantities 
directly related to the penetration distance of the beam into the gas, are best written 
down as a function of the following variable: 


=t z= C (5.1) 


has the dimension of the time and is a fixed label for a given particle within 
the pulse, provided this particle does not move relative to the beam head. For the 
beam head = 0, and for the tail-end, is simply the beam pulse duration. (See 
Figure 5.) 


The circuit equation has to be derived from the full set of Maxwell ’s equations. 
However, provided that a similarity assumption is made on the radial profile of the 
various charge and current distributions, the circuit equation takes a remarkably 
simple form. First, it has to be noticed that because the plasma electrons are mostly 
generated directly by the beam, ne, and thus , will have radial distribution close 
to that of I}. Second, as a consequence, the plasma current I, will also have a 
distribution similar to I, and fn , just like fe, will be independent of r. The 
natural radial distribution function to take (for all quantities) is that of a Bennett 
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profile (4.16), and in that case the circuit equation takes the form [34]: 


@ e? 
k= ay) ng CE) eng CH): 62 


The second term on the right is the current associated with the charge neutralization 
process and is thus important only at the beam head and neck. The third term is 
the plasma return current Ip. The last term is the displacement current which can 
be generally neglected, except at the very head of the beam. In this equation, e 
and m are calculated from the on-axis values of 


The exact form of the equation giving the plasma conductivity depends con- 
siderably upon the chemical composition and the pressure of the gas. The general 
form of the plasma electron density equation is as follow [35]: 

Sns c+ ate ate rf (5.3) 
The first term is the direct beam ionization rate where S ( ) is the Bethe’s stopping 
power, and w = 33:7 eV for air, the energy required to create one electron-ion 
pair. The second source term corresponds to ionization avalanche in the induced 
electrical field at the beam head. The third corresponds to conductivity electron 
loss by attachement to various molecules, and the last one to losses by molecular 
dissociative recombination. The rates ;, a and the coefficient , depend in 
general on the gas pressure, the electric field and the plasma temperature. In fact, 
for a real gas like air, a complete description would require a set of rate equations 
coupling the various populations of molecular and atomic nitrogen, oxygen, water, 
etc., in various states of ionization to the beam source terms. This leads to the 
development of complicated phenomenological models [86]. 


In order to calculate the conductivity , the collision frequency of the con- 
ductivity electrons with the plasma is also needed. This is the sum of two terms: 


= e Te7ng) + i (T.e7Ng) (5.4) 


The first term is the electron-neutral momentum transfer collision frequency, and 
the second one the electron-electron collision frequency [35]. For dilute plasmas, 
in which ng is less than a few per cent of ng, the molecular gas density, the second 
term can generally be neglected. 


Finally, a last equation giving the plasma temperature Te is also required. 
This is generally done by phenomenologically relating Te to E, B6] B7], or by 
constructing an equation relating the conductivity to Joule heating, thermal 
cooling by conduction and energy loss as free energy [B8]. 
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In conclusion, the plasma chemistry involved in the calculation of the conduc- 
tivity is generally rather complicated. However, for short pulses in dense gases 
such as air in the pressure range of 0.1 to 1 atmosphere, a good first approximation 
is obtained by assuming an average plasma temperature of Te | 2eV and keeping 
dissociative recombination with oxygen as the main cause of electron loss. In that 
case 

10° cm*s* and 18 ifn,s* (5.5) 


i.e., 486 Inst for ng 34 18cm?®. 


The four equations (4.18), (4.20), and constitute a full set of equa- 
tions providing a complete 1-dimensional description of the propagation of a beam 
pulse in a gas or plasma. The main assumptions leading to these equations are the 
paraxial approximation, the description of the plasma by a scalar conductivity and 
the similarity of all radial distributions. Whereas these coupled equations cannot 
be solved in general without a computer, their main features and implications for 
endo-atmospheric beam weapons can be derived by making some approximations. 


5.2 Conductivity generation and critical beam cur- 
rent 


When a high current beam pulse is injected into air at full atmospheric pressure, 
the avalanche ionization term can be neglected and (5.3) can be solved explicitly to 
give neas a functionof . With the further approximations (5.5), one can calculate 


m- 
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with r, the classical radius of the electron and e the electron charge. I, the 
conductivity generation current, which controls the rate of conductivity generation, 
is a critical current independent of pressure, and of beam parameters for 1, 
i.e., relativistic beams for which Bethe’s stopping power S ( ) is nearly constant. 
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The value of I, is typically of the order of 10 kA. {Į , the asymptotic value of 
m, increases with the beam radius and current, as well as with a reduction in 

atmospheric pressure. For pulses that are long relative to m Is =g, m will thus 

approach y , whereas for short pulses, or at the beam head reduces to 


Ip 


m ( ) 
ic 


(5.9) 


We will now make the drastic approximation that both a° and a® can be 
neglected in (4.18). This approximation is only valid in the quasistatic limit leading 
to Bennett’s or Nordsieck’s equations. However, this hypothesis corresponds 
approximately to the the situation in which a beam head leaves the exit window of 
an accelerator to penetrate into the atmosphere. Finally, we will assume a beam 
current I, with an infinitely short rise-time. Under these conditions, neglecting 
the displacement current, the circuit equation can be solved. With (5.9), the 
effective current is then 


1 a he 
p 
where 7 
= 2L— (5.11) 
Is 
and r 
a Ig 
= — 2— : 5.12 
a (5.12) 


It is now possible to calculate various quantities, and in particular the plasma 
current Ip which turns out to be maximum at t= p. The region where Ip is 
maximum can be defined as the beam neck. However, as p is calculated for a 
beam with infinitely short rise-time, it does not give the absolute position of the 
beam neck. Nevertheless, , can be used to calculate relative quantities, and in 
particular the neutralization fractions at the point where the beam is pinching. 
This is because in the beam head the plasma phenomena are dominant an to first 
approximation independent of the rise-time: the conductivity becomes smaller 
as the rise-time is reduced, but the inductive electric fields becomes greater, thus 
resulting in the same plasma current when I; reaches its maximum [39]. 


For = ,, the charge neutralization is 


1 1= 


fe = 1 1+ (5.13) 


At the pinch point, f, is thus equalto 1 1=e 0:63forl 1 and about one for 
1. The corresponding result for the current neutralization fraction is 


1 1= 


fp = = 14 (5.14) 


The maximum current neutralization fraction is thus small as long 1, but 
close to one for 1. Consequently, in order for a beam to be well pinched, 
which requires fm 1, one has to have < 1, or, explicitly 


1 
I, — 1 
B < a (5.15) 


AS Ig 10 kA, this is a rather strong limitation on the beam current. Indeed, for 
particle beam weapons, current in excess of I, are required for endo-atmospheric 
systems, and this creates some problems. 


These analytical results are in good agreement with detailed computer calcu- 
lations which do not assume a constant beam radius [40]. In these calculations the 
radius of the beam at the neck is well approximated by 


I 
a ie (5.16) 


which shows that for I; < Te, the neck radius is proportional to the beam current. 


The pinch condition is not the only one leading towards the requirement of 
small current neutralization fractions. We will see, for instance, that a beam is 
also more likely to be unstable when fw > 0:5. Therefore, it turns out that I, 
corresponds to a critical beam current setting an upper limit to the current of a 
beam pulse sent into the atmosphere. 


It is interesting to notice that the maximum current neutralization fraction cal- 
culated with is in reasonable agreement with the result of detailed computer 
calculations 40]. For example, in the case of 10 kA, 50 MeV beam pulse 
such as will be produced by the ATA accelerator which is under construction at 
the Lawrence Livermore Laboratory, we find fw p = 0:50, whereas the computer 
calculation predicts fwp = 0:4. In that case a= 05cm and b 10 cm; 
therefore L = 15and = 42. 


The maximum current neutralization fraction (5.14) is independent of pressure. 
But, as pressure is reduced, the ionization by avalanche can no longer be neglected. 
The result is then a lower bound to the maximum current neutralization 
fraction. For example, for the above mentioned ATA accelerator beam injected 
into air at 0.1 atmospheric pressure, the computer calculation give fnp = 08. 


The maximum value of the longitudinal electric field which appears during 
the self-pinching of the beam can be calculated from the circuit equation. In first 
approximation one finds: 


Ezmax = 7—g2e IA (5.17) 


In the head region where E , is maximum, ais of the order of b Thus: 


I 
E zmax = 2500- [V/m] (5.18) 


For beam currents of 10 to 100 kA and head radii of 3 to 30 cm, we find fields 
of the order of 10 MV/m. This maximum electric field is quite large and is 
approximately independent of the beam current. This is because the beam head 
radius becomes proportional to the beam current when Iz > I,, as it is shown by 
computer calculations [40]. 


In the body of a beam pulse long enough for m to reach its asymptotic value 
(5.7), the circuit equation becomes a linear equation with constant coefficients. It 
can then be solved easily. In particular, when neglecting te, the plasma current is 


I= exp Jfnp € (5.19) 


m 
For the previous examples, the plasma current decay time L ,, would be of 10 ns 
and 35 ns at the respective pressures of 1 and 0.1 atmosphere. Thus, if fa is large 
at the pinch point, the current neutralization fraction will remain almost constant 
throughout the pulse, provided its length is less than L m. 


A last case of importance is that of a beam pulse injected into a plasma of such 
a density that its conductivity remains constant. (This is in particular the case when 
e> min so that becomes the Spitzer conducivity [41].) In that situation, 
because the plasma density is high, plasma currents can start flowing right from 
the very beginning of the pulse. The current neutralization is therefore almost 
complete at the beam head, and with fap = 1 will give the plasma current 
throughout the pulse. Because of this large current neutralization, the beam head 
will not be pinched even though charge neutralization is complete. Therefore, a 
pulse injected into a plasma or preionized background will also have an expanded 
head, and beam pinching then occurs as the longitudinal plasma current decays. 


5.3 Beam head erosion and ohmic losses 


In all situations, whether the beam is injected into a neutral gas or plasma, the 
head of a pulse is not pinched and is followed by a neck in which there is a deep 
spike in the longitudinal electric field. The beam particles at the front of the pulse 
will thus spread out radially and get lost, and those in the neck will lose energy 
because of the longitudinal electric field. As their energy decreases, the particles 
in the neck will also experience more scattering, and both effects will combine 
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Head Neck Body 


Figure 5.3: Beam neck profile. The variations of the effective beam current I, , 
plasma return current Ip, kinetic energy K , longitudinal electric field E,, and 
beam radius a, are shown in the neck region which separates the well-pinched 
body from the freely-expanding head of the beam pulse. 


to increase the pinch radius. Consequently, the beam head expends continuously, 
and the neck region recesses progressively into the body of the pulse. (See Figure 
6.) 


This beam head erosion process has been studied both theoretically and ex- 
perimentally (34) {40}. These studies show that after an initial transient period, the 
erosion process results in an almost constant rate of decrease in the length of the 
pulse as a function of the propagation distance. 


The erosion rate can be estimated in a simple way by assuming that the particles 
in the neck will get lost when their kinetic energy has dropped from its initial value 
K down to some critical value K because of the ohmic losses in the longitudinal 
electric field. From and (4.20p, the erosion rate is thus 

d q d 


— (LI 5.20 
a (aK ve) (5.20) 


which can be integrated to give the propagation range z corresponding to the 
decrease tin the pulse duration: 


K K 
= —— .21 


Assuming K = K =2 as suggested by experiments [34], this gives 


A Ere (5.22) 


For a beam with K = 500MeV, I, = 10kA, and L = 1atthe neck, a beam pulse 
length loss of t= 10ns corresponds to a range of 600 m. Beam head erosion is 
thus a severe limitation to the propagation of short pulses. 


Moreover, the estimate has been obtained under the assumption that the 
particles in the beam neck do not move axially relative to the particles in the body 
of the pulse. This is a good approximation for electrons because in their case 

*  1at high energy. For protons, however, the condition * 1 requires very 
large energies, of the order of several GeV, in order to be satisfied. 


The ohmic losses are maximum in the head. However, for short pulses, the 
ohmic losses in the remainder of the pulse may also be important. From (4.7) 
and one can calculate the range of a particle within a pulse subject to these 
losses: 


= ->Š a -42 (5.23) 


The maximum range is of course obtained for particles in the tail for the pulse 
where m is maximum and £, minimum. 


The energy given up by ohmic losses primarily heats the plasma electrons. 
After the passage of the beam pulse, this energy contributes in forming a channel 
of hot air along the beam path. On the other hand, the energy lost in the form 
of radiations or particles scattered out of the beam do not contribute substantially 
into heating the air. The total energy deposited per unit length along the beam 
trajectory can be calculated by integrating the collision and ohmic terms in (4.20) 
over the length of the beam pulse: 

W ar 


L 
= ŽŽ Sc t+ 4L — (5.24) 
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where L and I, are evaluated at the end of the pulse. This result shows that, for a 
high current beam, the energy depositions by ohmic losses can be larger than the 
energy deposition by collisions. This happens when the gas pressure is low enough 
for collisions to be negligible, or when the pulse is short. Energy deposition can 
thus be enhanced by chopping a given pulse into a series of smaller pulses. If the 
individual pulses of the pulse train propagate independently from another, they 
will of course all be subject to erosion. The most efficient configuration for energy 
deposition will thus be a compromise. 
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Chapter 6 


Stability of propagating high-power 
beams 


6.1 General considerations on beam stability 


A crucial question concerning beam propagation is that of stability. Because of 
the substantial source of free energy represented by the kinetic energy of the beam, 
a variety of instabilities could be excited and amplified during propagation. For 
the purpose of using high energy particle beams as weapons, the problem is to 
find a set of beam parameters (energy, current, radius, pulse length, emittance, and 
energy spread) such that the beam can reach the target without being destroyed by 
the possible instabilities. This is rather difficult, especially because the existing 
theoretical models predict stable beam propagation in such a narrow range of 
conflicting parameters that only the actual testing of a beam will confirm whether 
these predictions were correct or not. 


The stability of beam-plasma systems is investigated by the standard perturba- 
tion method. If the initial perturbation of a stationary state of the system increases 
with time, the state is unstable under a perturbation of this type. Usually one seeks 
a solution of the system in the form 


fwo= f e)expilt K x) (6.1) 


where f is the deviation of any physical quantity from its stationary value. The 
relation between the complex frequency ! and the complex wave number X is by 
definition the dispersion relation 


D (!;K)= 0: (6.2) 
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A wave is said to be unstable, if for some real wave number k, a complex ! with 
a negative imaginary part is obtained from the dispersion relations, signifying 
growth in time of a spatially periodic disturbance: 


Im k) = 0 E = 
in (1) < 0 =) Instability (6.3) 
The absolute value of the imaginary part of ! for an unstable wave = 
Im (! ) is called the increment or growth rate (or its inverse the e-fold time) of 
the perturbation, because 
f(t) / exp t: (6.4) 


Instabilities in a beam-plasma system are primarily the result of the interaction 
between the beam and the plasma. The main parameters characterizing the plasma 
are its so-called plasma frequency (or Langmuir frequency) !,and__ the collision 
frequency of the plasma electrons. The plasma frequency is 

1g 
i = —— ne (6.5) 
oMe 
where ne is the electron number density and m e the electron mass. The plasma 
dielectric constant is then 


=i 2 (6.6) 


Similarly, as a plasma on its own, the beam is also characterized by its beam 
plasma-frequency (or beam Langmuir-frequency) 


= —— n (6.7) 


where n, is the beam particle number density, m the beam particles rest mass and 
their Lorentz factor. The on-axis beam-plasma frequency is related to the beam 
scale radius by (4.2), and therefore 


= 4—~— (6.8) 


Finally, with the definition of the plasma frequency, the magnetic diffusion time 
(4.8) can be written 

12 ¢ 
2E 
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(6.9) 
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and the charge neutralization time (4.3) 


e= > (6.10) 


There are numerous kinds of instabilities and many methods of classifying 
them. A rough phenomenological way is to divide them into two categories: 
macroscopic and microscopic. The next step is then to distinguish between var- 
ious characteristics, such as the fate of the perturbation (absolute or convective), 
the driving forces (electrostatic, electromagnetic, etc.), the relative importance 
of collisions (collisional/non-collisional, resistive/ non-resistive), the relative im- 
portance of spread in the velocity distributions or temperatures (hydrodynamic, 
kinetic), etc. There is no sharp distinction between different possible categories. 
In this paragraph we will simply introduce the main concepts involved. 


Macroscopic instabilities influence the spatial distribution of the beam. They 
are usually classified according to the geometry of the distorsion. For beams there 
are four main categories: (i) the sausage mode in which the beam contracts ar 
regular intervals; (ii) the hose or kink mode in which the beam oscillate sideways; 
(iii) the filamentation mode in which the beam breaks up into several filaments; 
and (iv) the ripple mode in which the beam is distorted by small-scale ripples on 
the surface. Macroinstabilities concern also other macroscopic degrees of freedom 
such as densities, hydrodynamic velocities, etc. They are connected with the flow- 
out of a plasma as a whole from one region into another. Macroscopic instabilities 
are also called magnetohydrodynamic, hydrodynamic, or simply low frequency 
instabilities. 


Microscopic instabilities do not necessarily induce a macroscopic motion of 
the plasma as a whole, but they can excite local fluctuations of density and electro- 
magnetic fields in the plasma. These velocity-space modes which do not appear to 
have direct effect on the beam will principally appear as an extra form of energy 
loss by which the energy of the beam can be transformed into powerful electro- 
magnetic radiations. The most important kind of microinstability is the so-called 
two-stream instability. This mode is found when one is studying the propagation 
of electromagnetic waves in a system consisting of two interpenetrating streams 
of particles. 


Some kinds of microinstabilities are directly connected to macroinstabilities. 
This is the case, for example, of filamentation which can be the macroscopic 
stage of a purely growing transverse electrostatic microinstability (the Weibel 
instability). In general, the growth of microinstabilities may lead to the onset of 
macroinstabilities, which may eventually destroy the beam. 


In a given reference frame, two types of instabilities can be distinguished 
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physically: convective instabilities, and absolute or non-convective instabilities 
[42]. Briefly stated, the essential point is the distinction as to whether an initially 
localized disturbance (on an infinite system) grows exponentially with time lo- 
cally (absolute instability) or ultimately decays because of the propagation of the 
growing disturbance away from the point of origin (convective instability). For a 
beam propagating along the z axis, the group velocity of an unstable perturbation 
with longitudinal wave number k, is given by 


@ 
Vg = ao : (6.11) 
If vy = v the beam velocity, the perturbation is absolute in the beam frame; if 
Vg = 0, the perturbation is absolute in the plasma; and if v, < v, the perturbation is 
convective. Clearly, for beam propagation, the worse instabilities are the absolute 
ones in the beam frame. Therefore, if an instability cannot be avoided, it should 
at least be convective in order not to completely hinder propagation. 


The effect on stability of beam particle collisions with the plasma can in general 
be neglected except when, as with bremsstrahlung losses of high energy electron 
beams, they contribute in broadening the beam velocity distribution. 


The collisions of the plasma electrons with the plasma molecules and ions often 
have an important effect on stability, and they can either increase or decrease the 
growthrates. When collisions are negligible, the instabilities are in the collisionless 
regime ( ! 0). In the collisional regime several distinctions can be made, in 
particular depending upon the relative values of the magnetic diffusion time m 
and the collision frequency [43]. For high-current beams propagating in air, 
one has m < 1 for pressures below a few Torr, and the instabilities are called 
non-resistive. For pressure above afew Torr , > 1, and the instabilities are of 
the resistive kind. In the former case, the key parameter is the dielectric constant 
and various time scale are possible for the growth rate, in the resistive case the 
time-scale is m. When e> ,,, the acting forces are predominantly electrostatic, 
whereas, in the resistive domain where m > .«, they are primarily magnetic. 


In describing the electromagnetic oscillations associated with instabilities, 
especially in the case of microinstabilities, one often uses the standard conventions 
used for electromagnetic wave propagation in wave guides or plasmas. If Æo and 
Bo are the unperturbed fields, and E, and B4, their respective perturbations, one 
uses the following terminology: 


B,= 0 electrostatic wave, 


B,€0 electromagnetic wave, 
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K kBo parallel wave, 


K? Bo perpendicular wave, 
K kE: longitudinal wave, 
K? Fy transverse wave. 


The hydrodynamic description of a beam-plasma system is strictly valid only 
in the limit of a monoenergetic beam penetrating a cold plasma. The velocity 
distribution of the beam (also called the beam temperature), and/or the velocity 
distributions of the plasma (the plasma ion and electron temperatures) can be taken 
into account by kinetic plasma theory. Neglecting the plasma temperature effects, 
the hydrodynamic regime of instability, as opposed to the kinetic regime, is defined 
by the following condition 


Xp vj> v hydrodynamic regime 


Xp vjýj< v kinetic regime (6.12) 


where v, = ! =k is the phase velocity of the wave, v the beam velocity and vthe 
RMS beam velocity spread. Introducing the Doppler shifted wave frequency 
=! K v; (6.13) 
the kinetic regime is then defined as 
j j< K y: (6.14) 
For self-pinched beams, the velocity distribution is the result of both the 


spread in beam energy ( W ) and in beam direction (emittance). The transverse 
and longitudinal components of the resulting velocity spread v are [44 


r 
V? Ig 
T= =~; 6.15 
Vk 4 2 1 W 
Isep y 6.16 
7 2 S (6.16) 


The combination of the effects of collisions with those of velocity spreads 
can considerably diminish the growth rate of instabilities. However, the most 
serious instabilities cannot completely be suppressed. Their effect will thus have 
to be minimized by the proper choice of beam parameters like current, shape, 
pulse length, deliberate velocity spreads, etc. In the following sections the major 
instabilities will be discussed together with possibles remedies. 
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6.2 Microinstabilities 


When a beam passes at high velocity through a plasma, its coupling to the electrons 
or ions can excite unstable electromagnetic oscillations. For high energy beams 
propagating in ionized air, the coupling to the ions is in general negligible relative 
to to the coupling to the electrons. The resulting microinstabilities are of several 
different kinds depending upon the type of electromagnetic wave excited. The 
two most important ones are the so-called two-stream instability, which refers 
to transverse or longitudinal electrostatic waves, and the Weibel instability (or 
filamentation) which is a transverse electromagnetic mode. 


The complete analysis of the various types of streaming instabilities is rather 
complicated because the beam velocity distributions, the plasma temperatures, 
the beam’s own magnetic field, etc, all have an effect on them. The theoretical 
analysis by analytical models is thus restricted to the most simple cases for which 
a linear perturbation approach is possible. The full analysis, including non-linear 
effects, requires computer simulations and experiments will be needed to check 
the validity of the codes [22] 23]. The following discussion will present the main 
known results for the above-mentioned instabilities in the case of high energy 
beams launched into the air in the pressure range of interest for particle beam 
weapons. 


6.2.1 Two-stream instability 


The dispersion relation of the two-stream instability of a monoenergetic beam 
interacting with a cool plasma, neglecting he beam’s own magnetic field, is [45] 


+ cs )=0: (6.17) 


where is the angle between the wave vector and the beam velocity 
K w=kvos = Ww (6.18) 


The unstable waves which are solutions to this dispersion relation are obtained by 
resolving it for ! and finding he roots which have a negative imaginary part. The 
peak growth occurs atk,v= !,andfor > !, itis given by 


r 
l l 


=In(!)=— (2s? +o2 ): (6.19) 


At relativistic energies, the fastest growing oscillations are those propagating 
almost perpendicular to the beam (sin 1). Such an angular dependence is 
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due to the fact that when 1, the transverse mass mg of the beam particles 
is much smaller than the longitudinal mass m °, and the oscillations which are 
transverse to the beam are easier to build up. The worst case will thus correspond 
to transverse waves for which 


T 
! z sh : (6.20) 


The most difficult mode to suppress corresponds to k} = 1=a because the beam 
has finite radius a. Therefore 


sin k =k; = (6.21) 
For air at atmospheric pressure, typical values for the plasma parameters of a 
10 kA, 1 GeV, 0.5 cm radius electron beam are: !;= 5 16°,!,= 6 16%, 
and = 47 16%. Thus sin = 0:01 and the efold distance corresponding to 
the maximum growth (6.20) is v= = 1m, showing that two-stream instability 
ia a major obstacle to the propagation of monoenergetic beams over significant 
distances. However, taking the self-fields of the beam into account, the two-stream 
instability can be totally suppressed, provided the parameters fall into the kinetic 
regime. The kinetic theory dispersion relation replacing in the case of 
transverse waves takes the form [42] 


ie d ew) 
—-_ e 0 : 6.22 
IC+i) ©? U k ae 
The integral over the beam velocity distribution f(v) has a singularity for ! = kv, 


which was first studied by Landau for the collisionless damping of electromagnetic 
waves in warm plasmas (Landau damping effect). In the case of the two-stream 
instability, a similar effect is associated with the beam velocity distribution and re- 
sults in the suppression of the instability [45]. This suppression becomes effective 
in the kinetic domain, which for two-stream instabilities sets in when 


jj mM(!)<k y : (6.23) 


The transverse velocity spread due to particle oscillations in the pinch field is 

given by (6.15). The stability boundary, which is also the worst-case kinetic 

regime growth, is from (6.20) and (6.21) 
r 


I> c I 
i, E a 6.24 
‘3 Ss rae 20) 

and [46] 

(6.25) 


The simple stability criterion predicts that the two-stream instability will be 
suppressed for high-current high energy beams propagating in air at pressures 
above a few Torr. This has been verified experimentally and has been 
confirmed by extensive computer calculations for beam energies up to 1 GeV and 
currents of the order of 10-100 kA [49]. 


Below the critical pressure related to the collision frequency bound (6.25), 
the two-stream instability is not suppressed, but its growth rate is considerably 
reduced by kinetic effects. The growth rates of the longitudinal and transverse 
electrostatic waves, respectively, are then 


i ! 
=> ; (6.26) 
x 2 i Vk 
and 
Fs v 
p V? 


and can be considerably lower than (6.20). 


6.2.2 Weibel (or micro-filamentation) instability 


The electromagnetic modes with K normal to v, Æ nearly parallel to v, and B 
normal to both v and F is a variety of the Weibel instability. This instability grows 
fastest in the absence of external or self-generated magnetic fields and is thus 
maximum for fm 1. It is an absolute instability which grows at perturbation 
centers of enhanced beam density which magnetically attract nearby beam particles 
and repel plasma electrons. Thus, the beam ultimately splits into filaments, each 
of which self-pinches. The computer simulations of this process show that 
the filamentation stage in which the effect of this microinstability reaches the 
macroscopic level, eventually ends up with the beam breaking up into separate 
filaments, which may recombine into a single dense beam, from which the return 
current is expelled. 


For a charge- and current-neutralized beam, in which the self-fields are nearly 
suppressed, the dispersion relation of the Weibel instability for low frequencies 


(|< !,) is BI] 
12 
1? = ae (6.28) 
1+ oe 


where k, is the transverse wave number. This is a purely growing mode and 
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in the case of a collisional plasma it becomes 


2 
13 = 4 228 Kc (6.29) 
Pp 


which shows that collisions enhance this kind of instability. In the kinetic regime, 
the growth of this instability is reduced, but kinetic effects alone are not sufficient 
to suppress it. This is because for such a transverse instability j j< j! jand there 
is no such effect as Landau damping in that case. The kinetic regime boundary 
gives the maximum growth rate of this instability. From (6.29) 


(6.30) 


From (6.8) and (6.15) the maximum e-fold time in the kinetic regime is finally 
1 1 
ous, G (6.31) 
(When £, = 1 the kinetic effects are absent and (6.29) applies.) The efold time 
is very short and would, for this absolute instability, destroy the beam in 
few nanoseconds. However, this is only the case for fm 1 and a presence of a 
sufficient non-zero magnetic field results in the stabilization of this mode. 


For collisionless plasmas, in the absence of an external magnetic field, the 
stability criterion is [53] 


2a G) UA 8 >z: (6.32) 
This condition indicates that for a charge neutral beam (fẹ = 1), a current neutral- 
ization fraction of less than 0.5 provides enough self-magnetic field to stabilize 
the micro-filamentation mode. For the case of collisional plasmas, there is no 
published analysis on the full effect of the self-fields. However, the criterion 
(6.32) is likely to still be applicable. This is, as is confirmed by recent theoretical 
work [54], because this mode is always stable when the magnetic field B gives 
a beam particle cyclotron frequency .= QB= m larger than !;,, and that the 


self-magnetic field is equivalent to = 2 !» (53). 


6.2.3 Discussion of microinstabilities 


Concerning the stability of a beam with respect to microinstabilies, the situation 
can now be summarized as follows: 
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At very low gas pressure, the plasma generated by the beam can be so 
week that the beam is only partially charge neutralized (fe < 1), and the 
conductivity so low that fm 0. In that case, when 1 > £ > 1= 7, 
the beam is pinched but the plasma electron density may be insufficient 
for microinstabilities to develop. With electron beams, the two-stream 
instability is avoided as long as the product of gas pressure and pulse length 
is low enough to keep secondary electrons from accumulating inside the 
beam |55]. This is sometimes referred to as the low pressure propagation 
window. For air, it corresponds to about 10 © Torr fora 10 kA beam [21] [55]. 
For such a beam, with an energy of 500 MeV, the range would be of the 
order of 200 km |55]. 


At intermediate pressures, the two-stream mode sets in, and the beam be- 
comes unstable. The range of the beams is then determined by the maximum 


growth rate given by (6.26) or (6.27). 


At pressures above a few Torr for beam currents in the 10-100 kA range, 
the two-stream instability is suppressed. This has been observed in many 
experiments — for a review see and for recent experiments [22] 47]. 
However, above the critical pressure, the rise in plasma conductivity enables 
a return current to flow and filamentation becomes possible as fa ! 1. 
Moreover, because of this increase in conductivity, other new instabilities, 
including macroinstabilities, also become possible. Therefore, the narrow 
stability window observed in air near 1 Torr is due to the suppression of 
microinstabilities by collisions and kinetic effects, and to the temporary 
absence of resistive macroinstabilities which set in at higher pressures. 


6.3 Macroinstabilities 


The small amplitude distorsion of the shape of a cylindrical beam can be described 
by giving the perturbed form of the beam surface: 


r=at rz; jt): (6.33) 


The standard method is to make a multipole expansion of the perturbation in a thin 
annulus as a function of the azimuthal angle 


X 
r= Im (z;gt)œosm : (6.34) 
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For small amplitude periodic perturbations, this can be Fourier analysed as a 
superposition of modes such that 
X 
r= A, woexspik z !t+m ): (6.35) 


m 


The mode m = Odisplays harmonic variations of beam radius with distance along 
the beam axis: this is sausage instability. The mode m = 1 represents transverse 
displacements of the beam cross-section without change in the form or in a beam 
characteristics other than the position of its center of mass: this is called sinuous, 
kink or hose instability. Higher values of m represent changes of cross-section 
from circular form: m = 2 gives an elliptic cross-section, m = 3a pyriform cross- 
section, etc. Modes withm > 1 are referred to as filamentation modes, because 
their growth leads towards the break-up of the beam into separate filaments. 


For a given azimuthal wave number m , various radial modes are possible. They 
are usually classified according to the degree 2n of polynomial eigenfunctions of 
the dispersion equation. For example, when m = 0, n = O corresponds to axial 
hollowing, n = 1 to standard sausaging, n = 2 to axial bunching, etc (7. A 
systematic classification of the modes excited by resistive macroinstabilities in a 
simplified helical orbit beam model has been given by Weinberg [58]. 


The primary concern with macroinstabilities is to determine their growth rate 
and the conditions for which they are sufficiently convective to insure that the 
growth will be of limited consequence for a sufficiently short beam pulse. The 
analysis is then simpler in the Doppler shifted frequency representation of the 
oscillations: 

expikz !t) = expi( zv i $! ) (6.36) 


In this representation, gives the position of a section of the beam pulse measured 
from the head of the beam (.1), and is the Doppler-shifted frequency (6.13). 


When is real, and ! complex, the dispersion relation D ( ;!) = 0 yields 
solutions giving the free growth of instabilities corresponding to initial value 
problems in the beam frame. These are potential instabilities developing during 
the flight of the beam towards its target. It is thus crucial that these instabilities 
are convective, i.e.,Re! ( ) 6 0. 


When ! is real, and complex, the same dispersion relation yields solutions 
to initial value problems in the accelerator frame. The convective nature of 
the instabilities or the limited duration for the acceleration period for absolute 
instabilities, will ensure that their effect will disappear after the pulse has left 
the accelerator. For practical reasons one usually studys macroinstabilities in 
laboratory experiments by deliberately introducing perturbations at the end of 
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the acceleration process [59]. In the following discussion, we will assume that 
the accelerator can be built in such a way that the pulse is not disturbed during 
acceleration. We will thus concentrate on instabilities affecting the flight of the 
beam towards its target. 


6.3.1 Electrostatic kink instability 


The analysis of the major macroinstabilities for propagation in air at pressure 
between 0.01 and 1 atmosphere, will show that only short pulses may propagate 
over sizable distances. Furthermore, because of the considerable difference in 
relativistic mass between transverse and longitudinal motions of the beam, insta- 
bilities involving transverse displacement wil grow fastest. As a result, the first 
order description of a beam pulse during the early stages of the growth of an 
instability, will not be that of a thin and flexible thread but of an axially rigid rod, 
the cross-section of which is undergoing distorsion. This is the so-called rigid 
beam approach. 


However, in the case of a beam propagating in a vacuum or a low density 
plasma, the thin thread approximation can be satisfactory. Let us consider, for 
example, the stability of a charged beam in a vacuum, and then the effects of 
a plasma on it 60]. In that case, the most serious instability is the m = 1 
electromagnetic kink mode in which the beam is distorted as a flat snake or a 
screw. The growth rate can be determined by analysing the lateral motion of a 
flexible cylinder of charged particles in a vacuum. The calculation is easiest in the 
beam frame, and after Lorentz transformation into the accelerator frame the result 


is s 
ka ka 
= ii— he) (6.37) 
Because of the finite radius of the beam, the growthrate = Im (!) = Im ( ) 
fork = k, real, has a maximum at ka = 0:6 , and a cut-off at ka= . Inthe 


rest frame of the beam, the effect of a plasma environment is to reduce the electric 
field by a factor equal to the dielectric constant : hence to reduce the growth 
rate by 1 >: When > kv, then becomes 


= ii? i kvne&): (6.38) 


This is the dispersion relation of the electrostatic kink instability in the hydrody- 
namic limit for a thin beam. It applies for beams propagating in collisional plasmas 
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under such conditions that the magnetic forces can be neglected. This is the case 
in air at pressures below 1 Torr. The maximum growth rate is at ka = 0:72 


1 
| E (6.39) 


In the kinetic regime, this instability is damped when IM A ( ) < k vx. This 
happens when 


m > =- i: (6.40) 


Let us consider now the stability problem of a rigid beam subject to transverse 
displacements. This is the m = 1mode and it is easiest to start by considering first 
a beam with a constant density profile and a sharp boundary, even though we will 
have to examine later the more realistic case of a beam with the Bennett profile. 
We further assume that the beam is fully charge neutralized by a collisional plasma 
of constant density with radius larger than the beams radius. 


If a beam is slightly displaced in the transverse direction, a surface charge 
density distribution with a cos azimuthal dependence will appear. Such a charge 
distribution produces an homogeneous dipolar electric field which results in an 
electrostatic restoring force directly proportional to the displacement. For a beam 
with a flat density profile, the equation of motion is simply 

2 
where y, is the position of the axis of the displaced, y, the position of the axis of 
the non-neutral plasma column which neutralizes the beam, and = = È + ve the 
total derivative. From the equation of continuity and Gauss’s law one can see that 
the plasma will move according to 


@ 1 


Yb = ZLE (Yo y) (6.41) 


a = 7 V“ yy: (6.42) 
Taking small perturbations of the form expikkz  !t), this couple of equation 
yields the dispersion relation 
2 
Uesg —@ (6.43) 
2°b 


This is the electrostatic kink instability. For real, ! is purely imaginary and 
the growth is absolute in the beam frame. Furthermore, the growthrate = i 
tends to infinity as ° approaches z 12. If we take for . the complete expression 


i! 


= (6.44) 


the dispersion relation can be put in the more familiar form [62] 


12 is 
Ie iy ‘ 
meeer A (6.45) 


Apart from the factor Z, it is identical to the dispersion relation for transverse (i.e., 
= 5) two-stream instability (6.17). In the collisionless limit ( = 0) this result 
is well known to particle accelerators specialists [61]. 


6.3.2 Electromagnetic kink (or “‘hose’’) instability 


In a low conductivity plasma, «is very large and reduces to stable os- 
cillations at frequency * = $!/. Physically, thy, this corresponds toa beam 
neutralized by infinitely heavy ions, so that the beam’s motion is simple harmonic 
oscillations about a cylinder of neutralizing charges which cannot move. This is 
obvious from which shows that when = 1 , y, is constant. At the other 
limit of high conductivity e becomes very small. In that case indicates 
that yp = yp, and this corresponds to the fact that when conductivity is very high, 
the plasma neutralizing the beam can follow it exactly, and there are no electro- 
static beam oscillations or instabilities. However, when the conductivity is high, 


magnetic forces have to be taken into account: 


When a beam is displaced transversely to its direction of propagation, the mo- 
tion of the magnetic field generated by the beam current Jg induces a longitudinal 
electric field. If v, = Ey is the velocity of the sideways beam displacement 


and B the azimuthal magnetic field, the Lorentz transformation of this magnetic 
field gives E, = v, B in the limit of small v,. This electric field generates a 
longitudinal plasma current J, = E,, which, if the plasma conductivity is high 
enough, generates an azimuthal magnetic field of sufficient strength to interfere 
with the beam current generated azimuthal magnetic field. (This plasma current 
should not be confused with the plasma return current, which we assume to be 
negligible for the moment, fm = 0.) As the displacement of the plasma does not 
coincide with that of the beam when the conductivity is finite, the total magnetic 
field resulting from Jz and J, will have an axis Ym different from the beam axis 
yb or the plasma current axis yp. This displacement of this magnetic field axis is 
directly related to the magnetic diffusion time. Thus 

a y: (6.46) 
When the conductivity is low (1.e., m small) the plasma induced by the lateral 
displacement is negligible and yn = y», the magnetic axis always corresponds 
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to the beam axis. On the contrary, when the conductivity is high (i.e., m large), 
or the plasma infinite in extent (i.e, L ! 1 ), Ym is constant and the magnetic 
field is “frozen.” When the magnetic field axis does not correspond with the beam 
axis, the particles in the beam are subjected to a restoring magnetic force which 
is equivalent to the force needed to drive the plasma current J,. This force can 
be determined from the effect on the beam of the dipolar magnetic field resulting 
from the differences in position between the beam and the magnetic field axis. 
This gives 

g — 1 2,2 . 

qn = 3 lp% y): (6.47) 
For periodic small amplitude oscillations, the system (6.46), gives the 
following dispersion relation [63] 64]. 


2 


UL, = r : (6.48) 


1 212 2 


2 “b 
This is the dispersion equation of the electromagnetic kink or hose instability of a 
rigid beam in a resistive plasma. It shows absolute instability in the beam frame 
and infinite growth for * ! + ?!2, This resistive hose instability has been 


extensively studied [57], [58], and to [69]. 


The dispersion relation (6.48) is similar to its electrostatic counter part (6.43), 
but corresponds to a different driving mechanism. In fact, for m < e the kink 
instability is of the electric kind and the dispersion relation is (6.43). On the other 
hand, in high conductivity plasmas m > eand the instability is of the magnetic 


kind and the dispersion relation is (6.48). In the limit of m = 0, the electric 
oscillations are stable with frequency * = $!/, and in the limit m = 1 the 
magnetic oscillations are stable with frequency *= 4 *!>. 


6.3.3 Macrostability of a beam penetrating a neutral gas 


In the case of a beam penetrating a neutral gas and generating it own plasma, the 
electric and magnetic kink modes are encountered successively as the conductivity 
rises from zero to a maximum. A model valid for arbitrary conductivity and 
combining the two instabilities is thus important. Such a model is obtained by 


combining (6.41), (6.42), (6.46), and (6.47): 
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Gave Z'o UM y+ i liy MK)=0; (6.49) 
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ev? T (Yp y)= 0; (6.50) 
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The resulting dispersion relation is [57] [69] 


(lL iLa)@ ile)= (6.52) 


As them = 1oscillations are stable in the limits = Oand = 1 ,the instability 
growth rate has a maximum at some finite value of the conductivity. The dispersion 
relation (6.52) shows that this maximum is obtained when 

aP — 

e= Lan b (6.53) 

c 
By comparison with and (5.12), one sees that this happens at a point very 
early in the beam pulse, and that m = 1 instabilities will have maximum growth 
in the neck region. 


In the head region of the beam, and m are direct functions of the beam 
current (5.9). Therefore, in order to minimize the length of the neck region over 
which the growth is largest, the conductivity generation by the beam should be as 
large as possible, requiring 

ip Te. 3 (6.54) 


However, this condition implies high plasma return currents: the hose instability 
should therefore be examined in the fa 6 Ocase as well. 


A non-zero plasma return current has essentially no effect on the electrostatic 
forces, and thus on the electric kink instability. On the other hand, the magnetic 
forces are directly affected, and the first effect of a return current is to diminish the 
magnetic force in by the factor (1 f ) which is equivalent to replacing 
the current I, used to calculate ! 2 with by Iz. But there is also an additional 
effect: as anti-parallel currents tend to repel each other, when the beam is displaced 
relative to the plasma channel in which the return current flows, the interaction 
between the beam current and the return current tends to increase the displacement. 
This is the so-called self-hose effect which results in a supplementary force on the 
beam proportional to Is Ip = fm If. Including those two effects, the equation of 
motion becomes 

g 1 242 1 242 
Y= 3 lE WIA f)t+ 5 “!symfa 
cd (6.55) 


= $7)? y0 ff) M 
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In (6.52) this results into the replacement of !2 by !7 (L {f ) in the denomi- 
nator on the right hand side [57], and for n > e 


= 
Le => g -~ (6.56) 


With the effect of the self-hose included, the growth rate of the magnetic hose 
instability is clearly worse. In particular contrary to (6.48), the growth rate at 

= 0 is now non-zero. This implies that when f, € 0, even a very slow 
transverse displacement may result in the disruption of the beam. Furthermore, as 
fm tends towards one, the growth rate increases and the position of the pole moves 
towards smaller frequencies. 


The electric and magnetic hose modes for rigid oscillations of flat profile 
cylindrical beams show very bad stability properties. The consequences of various 
damping effects on these instabilities will be discussed below. Before that, we 
will examine them 6 1 modes with the same flat profile rigid beam model, but 
only for resistive modes in the limitof m > & 


Form = 0, the sausage mode dispersion relation is as follows 


1 “ed 2+ * 
i=, = Se : (6.57) 
6 a2 R) 7? 
The main difference between this dispersion relation and the hose dispersion 
relation (6.56) is that the sausage mode is unstable in a limited frequency range 
only: for fn 0, between !, and 2!,. This will enable the kinetic effects to 
stabilize this instability provided that, as can be seen from (6.57), 
fa< 4: (6.58) 
Form > 1, there are numerous unstable modes possible. However, in the limit 
of small fa , the main modes obey the approximate dispertion relation [58] 
m? 1E “tet? ? 
i! m ~ma 2 : (6.59) 
a“ 
In the casem = 1 one recognizes the hose mode, and for m > 1, similarly to the 
sausage mode, instabilities exist only in a limited frequency range, which becomes 
narrower in proportion when m ° increases. 


The essential difference between the hose mode and the m 6 1 modes is 
that when m = 1a low frequency disturbance can be produced without internal 
compression or distorsion: only a simple transverse displacement is required 
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[64]. Because there is no change in internal pressure to produce a restoring 
force, hose instability appears at an arbitrary low (Doppler shifted) frequency. 
Furthermore, the stabilizing effects of possible spreads in the beam velocities 
vanish in that limit. For them 6 1 modes, instability potentially appears at finite 
frequency, but is strongly suppressed by kinetic effects. We will thus concentrate 
on these effects on the hose mode only, assuming that coping with this worst case 
mode will be sufficient. 


6.3.4 Macrostability of beams with rounded radial profiles 


In examining now the stability of more realistic beam plasma models we will see 
that the instabilities are somehow not so bad as for idealized models such as the 
sharp boundary, flat profile, cylindrical beam studied so far. We will first look at 
the effect of rounded beam profiles — specifically of the Bennett profile (4.16) — 
and then at the effect of spreads in velocity distribution, mainly arising from the 
particle’s oscillations in the pinch field. 


Let us rewrite the hose dispersion relation for fw = O in the following form 
2 
ie = -zz 2? (6.60) 
“b 
! is the betatron frequency, equation (4.15), which for a beam with a flat current 
profile it is simply 
! el? : (6.61) 
In the case of a beam with the Bennett profile (4.16), the betatron frequency is 
not constant but uniformly distributed between zero and the maximum ! m. The 
force equation has then to be averaged over the Bennett distribution and the 
resulting effective betatron frequency is [64] 


1 
=a aye (6.62) 
Similarly, for the magnetic axis diffusion equation (6.46), the dipole magnetic dif- 
fusion time must be recalculated by properly averaging over the beam and plasma 
conductivity profiles. Assuming the beam current and the plasma conductivity to 
Both have Bennett profiles with scale radius a for the beam and b for the plasma, 


one finds 


l 2 


Pi- ET A (6.63) 
Saa = ae i 
where = kæa. In the case where a = h this is 
1 
a> 8 m (6.64) 


52 


and by comparison with (6.63) one sees that a > = , when b> a, so that the 
beam is more stable in the case where the conductivity profile is broader than the 
beam profile [68]. 


The dispersion relation (6.60) has been derived in the hydrodynamic limit. 
However, from (6.60) and (6.15), 
a J j _ 
j j=! P < Km in v2 = ! : (6.65) 
Lae ol a 
Thus, according to (6.14), kinetic effects cannot be ignored for the hose insta- 
bility. The simple dispersion relation should therefore be replaced by the 


expression 
om 2 

ilg= El) ad (6.66) 

0 ! b 

where f (! *) is a suitable distribution function. In principle this distribution func- 
tion could be derived from the Vlasov kinetic plasma theory. But an exact solution 
to this problem has still not been found. However, two different phenomenological 
approaches [64 [65], which agree with computer simulations [68], give plausible 
results with the following function 


12 12 
et ! 
Oe? E (6.67) 
For fa = 0, the resulting dispersion relation is 
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In comparison with (6.60), this dispersion relation shows much less serious insta- 
bility problems: 


The growth rate has a cut-off at = !,, and there is no poleat = ! . 
Instead, for = 0:52! ,,, the growth rate has a maximum 
0:69 
max — : (6.69) 


d 


The instability, from absolute in the beam frame, becomes convective. 
Therefore, as a perturbation of the beam grows, it will at the same time 
move backwards into the pulse. For a perturbation of amplitude y (0) gen- 
erated at the head of the beam, a saddle point analysis of the dispersion 


53 


relation shows that the growth of the hose is such that at the tail of the 
beam pulse its amplitude will be 


t 
vi = y0) 1+ — : (6.70) 


The existence of a maximum and of a cut-off in have been verified 
experimentally [59]. Similarly, other experiments [47] have also demonstrated the 
convective nature of the hose instability. In these experiments, the hose instability 
leads to an erosion of the beam tail. This is because, for a beam of finite duration, 
the disturbance is maximum at the end of the pulse. 


The significance of is that for a pulse of finite duration, a disturbance 
ultimately disappears. Thus, if a beam pulse has a duration of the order of a, it 
will be able to propagate over large distances. In practice, the problem is that ą 
is only of a few nanoseconds for a high current beam injected into the atmosphere 
at normal pressures. Therefore, stable beam propagation is restricted to very short 
pulses. 


A further difficulty with the hose instability is that its maximum growth and 
its convective nature depend on the extend of current neutralization of the beam. 
It is found, for instance, that even a small amount of current neutralization (i.e., 
fn 0:1) is strongly destabilizing, particularly at low frequences, and leads to 
absolute instability as seen by the beam [65]. The hose instability in the presence 
of current neutralization has been studied in the frame work of he Vlasov theory 
[68]. The results show that for a current neutralization fraction larger than 0.5 the 
return current driven self-hose becomes the dominant destabilizing mechanism. 


The general properties of the hose instability can be transposed to the case of 
them 6 1 macroinstabilities and this has been verified in a number of cases. For 
example, it has been found that the sausage mode is stable provided that 


fn co , (6.71) 
if © 


Similarly, general considerations show that filamentation can be avoided if 


— (1 =1: 72 
a f) (6.72) 


6.3.5 Discussion of macroinstabilities 


In conclusion, the stability with respect to macroinstability seems to be insured 
for beams of pulse length of the order ot m provided that fw is small. In practical 
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cases, these requirements may be somewhat less stringent. For instance, the 
effect of various damping factors (such as variations of beam current during the 
pulse, energy spreads, smooth radial and longitudinal beam profiles, finite pulse 
length, etc [2I]) cannot be calculated easily. Experiments such as those that 
will be possible with the high-current high-energy accelerators which are under 
construction in Los Alamos [22], Livermore [23] and in the Soviet Union [72] are 
thus crucial. 


Finally, because of the need for short pulses, actual beam weapon systems 
will have to use trains of small pulses in order to send sufficient energy towards 
the target. This creates various problems, because pulses within a train will have 
to propagate in the plasma background generated by the previous pulses. The 
analysis of these problems, together with that of the boring of reduced density 
plasma channels, will be part of the conceptual design of practical particle beam 
weapons. 
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Chapter 7 


Discussion of theoretical prospect 


We have examined the main physical problems involved in the propagation of 
particle beams for possible exo-atmospheric or endo-atmospheric beam weapon 
systems. 


There is no problem of principle with the propagation of neutral particle 
beam in outer space. Systems with adequate characteristics are currently in the 
development stage [9]. For short range exo-atmospheric systems, i.e., of the 
order of a few hundreds of kilometers, charged beams or plasmoid beams might 
eventually be used. 


For land-based systems, the main problem arises because of the conflicting 
nature of the following set of conditions: 


1. The beam power has to be larger than the Nordsieck power (4.22) in order 
for the beam to propagate over sizable distances; 


2. the beam current should be less than the critical current (5.8), and the current 
neutralization fraction as small as possible (5.15), so that the beam is well 
pinched; 


3. the beam pulses have to be long enough so that they will not be completely 
eroded before reaching their target (5.22). Furthermore, the beam should 
be relativistic enough ( > 10) for the steady state erosion hypothesis to be 
valid; 


4. the pulse duration has to be of the order of the dipole magnetic diffusion 
time so that propagation is not disrupted by macroinstabilities (6.70); 
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5. in order to keep the growth of instabilities at a minimum, the current neu- 
tralization should be minimized (6.32), (6.56), (6.58), (6.71), and and 
the conductivity should rise fast enough to avoid excessive growth in the 


head (6.54); 


6. the beam radius, emittance, energy and current should match the Bennett 
pinch radius and should be compatible with the requirements of beam 
stability and hole boring. 


These conditions cannot easily be simultaneously satisfied. This is primarily 
because the pulse length is set by magnetic diffusion time. For a given beam 
radius, this time can only be increased by increasing the beam current, and thus 
producing an undesirable large current-neutralization fraction. 


However, within rather tight limits, a set of acceptable beam parameters seems 
to exist. In order to find out if these conditions are really satisfactory for beam 
weapon applications, careful experiments are needed. This concluding sentence 
was written in November 1982, together with the addition al remark: These 
crucial experiments will be done with accelerators currently under construction 
in the United States and probably in the USSR [72]. Within a couple 
of years, at most, the final answer should be known. What can we add to this 
conclusion, more then twenty years after it was written? 


First: That the laws of physics have not changed, and that the theoretical 
analysis done in previous chapters remains valid. In fact, having followed 
the literature over the past twenty years, it appears that nothing really new 
and important has been found and published. In particular, the most difficult 
theoretical problems addressed in this report — the questions related to of 
the stability of propagating a charged-particle beam through a pre-existing 
or beam-generated plasma — are still discussed in recent papers by referring 
to the seminal work done during the 1960s to the 1980s.' This is not to 
say that no progress has been made, but that the basic physical processes at 
work in high-intensity high-energy beam propagation in plasmas appear to 
have been properly recognized and analysed during these decades. 


Second: That a major step in the theoretical understanding of beam-propagation 
stability was made between the early analyses, such a those of Steven Wein- 
berg [58], and those of Edward P. Lee [64], which showed that oversimplified 
models were likely to predict far too pessimistic behaviors in comparison 


'See, for example, references number 20 to 25 cited in a recent paper where the resistive hose 
instability is reinvestigated for low-collisional plasmas {73}. 
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to more realistic ones. This means that an essential step after accelerating 
a beam pulse is that of “conditioning,” that is of shaping and smoothing its 
spatial an possibly temporal extent in such a way to avoid the onset and 
growth of instabilities during propagation.” 


Third: That technological advances are constantly changing the context in 
which specific systems have to be assessed, both from a theoretical then 
from an experimental point of view [74]. For example, the possibility of the 
use of antihydrogen beams in exoatmospheric systems came to the author’s 
mind in 1983 when a former CERN colleague told him about Los Alamos 
management’s interest to hire him for antimatter work at their laboratory 
[75]. This lead him to realize that antimatter was the sole portable source of 
muons, which are the main byproducts of the annihilation of antiprotons with 
matter [76]. This means that if muons from such a source could be “cooled”? 
and accelerated, the possibility of exploiting the high-range of muons beams 
(see Figures 4.2 and 4.3) would become theoretically possible. In fact, as 
will be seen in Section 10.5, the problem of generating and cooling intense 
beams of muons is the subject of vigorous research since a few years. While 
this effort may not succeed in producing muon-beams sufficiently intense to 
propagate in self-pinched mode over significant distances, the development 
of muon technology for fundamental research will help assessing whether 
muon beams, and beams of their related neutrino decay products, may have 
other practical applications [77] [78]. 


This is the plasma-physical counterpart of the hydrodynamical need to smooth boundaries 
and remove all unwanted discontinuities in order to avoid drag, turbulence, and instabilities in the 
motion of ships, aircrafts, or missiles. 

3“Cooling” means decreasing the kinetic energy dispersion of a bunch of particles. 
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Chapter 8 


Discussion of beam propagation 
experiments 


The discussion of the experimental prospect of particle beam technology for exoat- 
mospheric applications is quite different from that of endoatmospheric applica- 
tions.' Indeed: 


Exoatmospheric beam technology is very similar to the relatively low- 
intensity (mA), relatively high-energy (GeV), beam technology used in 
numerous, officially non-military, national and international nuclear and 
elementary-particle research laboratories. Moreover, the typical accelerator 
and beam steering technologies required for outer-space beam weapons is 
very similar to those of emerging applications of accelerators such as cancer 
therapy, breeding of fissile and fusion materials, transmutation of nuclear 
waste, breeding of antimatter [80], or proton-radiography of the detona- 
tion of conventional explosives and nuclear weapon’s primaries, and 
p.84—89], as well as to that of next-generation research accelerators 
such as the Large hadron Collider (LHC) at CERN and the future linear- 
collider to be built (most probably in the US) as a fully international effort 
[83]. This means that most components of space-based beam weapons can 
be developed and tested in civilian laboratories (Section 9.1), and that the 
role of military laboratories is to develop and test systems which integrate 
these components, and merge them together with the application specific 
technologies related to their use in an outer-space battle-field environment, 


'A summary of the early efforts on particle beam weapons research, from World War II to 
1980, as well as some comments on Soviet particle beam research, are given in a comprehensive 
fact-sheet published by the U.S. Department of Defense ‘ 
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rather than as research tools in a laboratory (Section 9.2). Finally, most 
experiments on exoatmospheric systems can be done on the ground, e.g., 
using evacuated pipes for propagating the beam, and without requiring safety 
precautions which are dramatically more stringent than those required by 
the operation of typical accelerators used for scientific research. 


Endoatmospheric beam technology is very similar to the relatively high- 
intensity (kA), relatively low-energy (MeV), beam technology used in mil- 
itary laboratories for the flash x-ray radiography of the detonation of con- 
ventional explosives and nuclear weapon’s primaries [82] p.81—84], the 
generation of high-power radio-frequency or free-electron-laser beams, and 
the simulation of nuclear weapons effects [29], with only few civilian appli- 
cations (such as industrial radiography, electron-beam welding, and possibly 
future inertial confinement fusion drivers and very high-energy particle ac- 
celerators [129]). Therefore, the development and demonstration of this 
technology is mostly done in military and national laboratories, as will be 
seen in sections 10.1 to 10.4, — even though there is significant overlap 
between these developments and those of conventional research accelera- 
tors, as is shown by the fact that both type of technologies are generally 
discussed at the same national and international scientific conferences, and 
published in the same technical journals. Finally, an additional reason for 
locating these experiments in military laboratories is that even relatively 
modest endoatmospheric systems create special hazards because the total 
energy in typical beam pulses can easily be equivalent to the energy content 
of gram- to kilogram-amounts of high-explosives. 


The purpose of this report is to review the physics of high-intensity high- 
energy particle beam propagation. The discussion of exoatmospheric systems 
could therefore be very brief since (as noticed in section 3.1) there is in principle 
no scientific obstacle to propagating a neutral beam over long distances in the 
near-vacuum of outer-space. We will nevertheless address some of the critical 
issues related to these systems, even though they are more of a technological 
than physical nature. This is because the development of exoatmospheric beam 
weapons is in many ways linked to those of endoatmospheric ones, and because 
their development illustrates the very strong interdependence which characterizes 
civilian and military particle accelerator research, as well as the importance of 
informal international collaborations and scientific exchanges (including between 
“enemies”’) in relation to the development of advanced weapons systems. 
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Chapter 9 


Neutral particle beam propagation 
experiments 


9.1 Neutral particle beam technology development 


According to section 3.1, what is required to focus a beam on a 1 m target at a 
distance of 1000 km is a beam with an initial radius of 20 cm and an emittance 
of 2. 10’m_ rad.As a charged particle is needed in order to be accelerated by 
electromagnetic means, the general technique is to start from a negative-ion source 
producing a low-energy, low-emittance beam of H or D atoms which is then 
injected into an accelerator whose main quality is to increase the energy of the 
ions to about 100 MeV without unduly increasing the emittance of the beam. At 
the exit of the accelerator this high-brightness beam is very precisely focused on a 
target by a magnetic beam-steering optics, and then passed through a neutralizing 
cell in order to remove the extra electron (without increasing the beam emittance) 
so that the final beam is neutral. 


Most of the basic technologies necessary for building such a system have been 
invented, built, and tested between the late 1970s and the early 1990s. This is 
what will be summarized in the next subsections: 

‘Chapter VII and VIII of the Los Alamos Accelerator Technology Division 1978 progress 


report is describing a number of features associated with a system based on a 20-100 mA, 50-500 
MeV beam with such an emittance p.24-77]. 
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9.1.1 The ion source. 


The first advance which made the concept of a neutral particle beam weapon 
scientifically feasible was the invention of a sufficient low-emittance negative 
ion source by the group of V.G. Dudnikov in the USSR [85]. This pulsed 
surface-plasma ion source was quickly adopted world-wide, and Los Alamos 
soon demonstrated its superior performances relative to conventional sources 


p.117], [87]. 


9.1.2 The injector. 


The second major advance which made the concept of a neutral particle beam 
weapon scientifically feasible was again a Soviet invention. It solved the problem 
of injecting into a linear accelerator a low-energy low-emittance beam while main- 
taining its brightness. It is so-called “radiofrequency quadrupole” (RFQ), which, 
in combination with the Dudnikov high-brightness H source, led to speculations 
in the years 1978-1980 that the USSR might have had a considerable advance in 
particle beam weapons technology. However, while development of RFQ tech- 
nology started in Russia in 1970 already, it only reached maturity around 1977, 
about the time when Los Alamos started its own RFQ programme [88]. 


In fact, it is in 1978 that a new Accelerator Technology Division (AT), headed 
by Edward A. Knapp, was formed at Los Alamos. The first step in the development 
of a neutral particle beam system for possible deployment in space, code-named 
“White Horse,’ was to build a 5 MeV accelerator test stand, which could be 
scaled up to about 50 MeV for outer-space re-entry-vehicle/decoy discrimination, 
and later to about 500 MeV for boost-phase or mid-course intercept of ballistic 
missiles. This accelerator test stand (ATS) was to integrate in a test-bed a Dudnikov 
H source, a RFQ to reach an energy of about 2 MeV, a drift-tube linear accelerator 
section to increase the beam energy to 5 MeV, and beam diagnostic equipments. 
This relatively small and low-cost system was to be supplemented by several 
collaborative undertaking on related military or civilian projects. This is why 
Los Alamos was glad to accept a proposal form CERN, reported in the 1980 AT 
division progress report p.36]: 


“(CERN) in Geneva, Switzerland, has asked the AT division to coop- 
erate in the design and construction of an RFQ linac. This linac would 
replace the [...] injector [...] used with the “old” linac for the CERN 
proton synchrotron. [...] The plan is to carry out the mechanical 
design and assemblies at CERN, and then send them to Los Alamos 
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for machining of the pole tips.” 


The reason why CERN needed Los Alamos for building an RFQ is that the 
precise machining of its pole tips required a computer-controlled milling-machine 
of a type that was only available in a nuclear-weapons laboratory such as Los 
Alamos. Reciprocally, Los Alamos needed CERN as it provided a test-bed for 
coupling an advanced RFQ injector to a high-energy linac, which, moreover, 
happened to be the most recent and advanced in the world at the time. 


Therefore, while the first successful high-power operation of a prototype RFQ 
with beam was obtained in February 1980, the first operational RFQ built at Los 
Alamos was shipped to CERN and installed on its new linac, while the second to 
become operational was installed on the Los Alamos ATS on November 16, 1982, 


[89] p.53]. 


9.1.3 The accelerator. 


A major achievement of the Los Alamos accelerator technology division has 
been to reach, in 1979, a quantitative understanding of emittance growth in linear 
accelerators [90], showing that there is a lower limit to the emittance of such accel- 
erators [91]. This understanding had been made possible by numerous exchanges 
between Los Alamos and many other institutions, including the CERN laboratory 
in Geneva. As an example, we cite a few sentences of the section on emittance 
growth of the 1978 Los Alamos accelerator technology division progress report 
(92) p.63,64,71]: 


“The new 50-MeV injector linac of CERN is the latest of this type 
proton linear accelerator to be commissioned in the world. [...] Pre- 
liminary results from CERN during the early summer 1978, using 
the first tank at 10 MeV, indicated larger than anticipated emittance 
growth. Knowing our urgent need for verification of the design ap- 
proach for the FMIT project, the CERN Linac Group graciously in- 
vited our participation in the commissioning of the full machine. An 
on-site collaboration afforded us the opportunity to reach a level of 
understanding impossible to achieve at conferences or by letter. [...] 
We very much appreciate the complete openness, candor, and hospi- 
tality of the whole CERN linac group. Their willingness to let us use 


We stress that these quotes are from a 47-page long chapter entirely dedicated to the discussion 
of the design of an accelerator and a beam-steering system to be used in a particle beam weapon. 
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their raw data, before they had time to publish their results, was the 
key that enabled us to answer our basic questions in the urgent time 
scale imposed by current projects.” 


9.1.4 The beam focusing and steering optics. 


A major issue in a space-based neutral beam system is the ability to focus the beam 
into an acceptably small spot at some distance, and to provide a steering capability 
for aiming the beam. Using the ray-tracing computer programme TURTLE (a 
typical particle beam optics development tool devised by a collaboration of ac- 
celerator physicists from FNAL, SLAC, and CERN), the Los Alamos accelerator 
group was able to show in 1978 that a 50-meter long system of bending magnets 
and solenoids was able to focus 91% of the beam within the desired radius, starting 
with an initial radius of 0.2 m and an emittance of 2 10’ m rad] p.52-62]. 
An optical system with similar characteristics was ultimately built and successfully 
tested in 1990 at the Argonne National Laboratory near Chicago. 


9.1.5 The neutralizing cell. 


The final part of a neutral particle beam accelerator system, the neutralizing cell 
in which the extra electron is detached from the accelerated ion to form a beam 
of neutral atoms, is possibly the most controversial component of the whole 
system. This is because simple techniques (such as passing the beam through 
a foil or a gas to strip-off the extra electron) are likely to significantly increase 
the emittance of the beam and therefore to unduly decrease its brightness, at 
least for the most demanding applications of outer-space weapon systems. The 
most promising technique is therefore that of laser photodetachment, which is 
potentially 100% efficient rather than limited to about 55% as with collisional 
processes. High efficiency in outer-space is very important since any degree 
of incomplete neutralization of the beam could result in charging the orbiting 
accelerator, ultimately causing arcs that could lead to self-destruction. Moreover, 
it could be that by means of a multiphoton detachment process the excess electron 
could be removed without increasing beam emittance [93]. However, the kind of 
laser required for this purpose would add a considerable level of complexity to 
a space-based particle beam device, so that the whole subject of neutralization is 
still an open issue, despite constant progress [PA [95]. 
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9.1.6 Summary 


The developments presented in this section clearly show that some kind of neutral 
particle beam weapons can be built. Therefore, they established the scientific and 
technical feasibility of such systems. However, they does not prove their practical 
engineering feasibility, namely they do not demonstrate that a complete system 
(including power supplies, cooling, and many ancillary systems) can be put in orbit 
and operated in outer-space. This is where the specific developments described in 
the next two sections come in. 


9.2 BEAR and GTA at Los Alamos National Labo- 
ratory 


Until 1983 the development of neutral particle beam weapon’s technology and 
systems was made in the form of small, relatively low-cost, and mostly “paper” 
research programs. This was because (as shown in the previous section) the much 
more costly civilian research programs related to fundamental research could be 
used to develop and build all the key components of such systems, and that even 
international laboratories such as CERN, located in a neutral non-nuclear-weapon 
state, could be used for testing key component in full-scale accelerators.° 


For this reason, the accelerator test stand (ATS), extensively referred to in the 
previous section, was sufficient to provide first hand experience on the most crucial 
components of a neutral particle beam system, and was even able to demonstrate 
reliable operation at 170 mA and 5 MeV, making it the brightest high-current H 
beam in the world [98] p.3]. 


Nevertheless, with the advent of President Reagan’s Strategic Defense Initia- 
tive, neutral particle beam weapons were to be evaluate in competition with other 
directed-energy systems such as high-energy lasers and nuclear-pumped x-ray 
lasers. This lead to the definition of a comprehensive, national neutral particle 
beam research and development program lead by the Los Alamos laboratory, de- 
scribed in a September-1986 brochure [98], and related technical reports [99]. 
The objectives of this program were: (1) to provide the necessary basis for a de- 
cision by 1992 to build a space-based neutral particle beam system to be used as 
a decoy/warhead discriminator (near term goal); (ii) to develop the technology 


3The same will happened in the mid-1980s, and still continue today, with the development 
of antimatter technology where U.S. scientists from weapons-laboratories have essentially free, 
unlimited access to all of CERN’s antimatter facilities [80][96] 07]. 
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in stages to ultimately build a neutral particle beam weapon capable to provide a 
rapid hard kill of enemy warheads (far term goal). To meet these objectives and 
more directly the first objective the essential elements of the program were: 


1. The Integrated Space Experiment I (ISE-1) and 
the Beam Experiment Aboard Rocket (BEAR); 


2. The Ground Test Accelerator (GTA) and 
the Technology Program. 


The near-term goal was therefore to use the space shuttle to launch, in 1991, a 
50-MeV-accelerator-based system into space (ISE-1), a quite ambitious objective. 
This was to be preceded by the suborbital launch, in 1987, of a 1-MeV accelerator 
by a single-stage Aries rocket. In this context, the GTA was to be representative of 
what was actually to be placed into orbit in ISE-1, while the BEAR payload was 
basically to be a ruggedized and less powerful version of the ATS, packed into a 
cylindrical volume of 1 m diameter and 7 meter length, together with diagnostic 
instrumentation. 


However, as a consequence of the January 28, 1986 Challenger space shuttle 
disaster and funding constraints, the Integrated Space Experiment (ISE-1) was 
cancelled, and Los Alamos was asked in December 1987 to complete BEAR on 
an accelerated schedule with limited funding p.25]. The BEAR payload was 
launched 200 km into space on July 13, 1989, and several accounts, e.g., reference 
[100], as well as the final report, i.e., reference [01], described the experiment as 
a success: 


“The US Department of Defense’s Strategic Defense Initiative Orga- 
nization is sponsoring the development of neutral particle beam (NPB) 
technology for strategic defense applications. The first step in demon- 
strating the functioning of an NPB in space was the development and 
launch of the Beam Experiments Aboard a Rocket (BEAR) in New 
Mexico in July 1989. A government, laboratory, and industrial team, 
under the technical coordination of Los Alamos National Laboratory, 
designed, developed, and tested the BEAR payload. The primary 
objective of BEAR was the operation of an NPB accelerator in space. 
The payload was also designed to study (1) the effects on the space ve- 
hicle of emitting an NPB and associated charged beams into the space 
environment; (2) the propagation and attenuation characteristics of 
an NPB in space; (3) the dynamics of the charged particle compo- 
nents of the beam in the geomagnetic field; (4) the effects of neutral 
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effluents from the vehicle; and (5) any anomalous or unanticipated 
phenomena associated with operating an NPB in the space environ- 
ment. The BEAR experiment successfully demonstrated operation of 
an NPB accelerator and propagation of the neutral beam as predicted 
in space, obtained first-of-a-kind NPB physics data, and demonstrated 
the ability of the BEAR accelerator to survive recovery and to continue 
operating normally. No unanticipated phenomena were encountered 
that would significantly delay further development of NPB technology 
for defensive, space-based weapon systems.” [101 


On the ground, despite the cancellation of the first (and thus subsequent) 
integrated space experiment(s), construction of the GTA proceeded more or less 
on schedule, and the accelerator was commissioned in 1992, producing a 24-MeV, 
50mA beam with a 2% duty factor [102]. Consistent with the design report, [99], 
and other reports, e.g., [98] [03], work on GTA and related equipments, such 
as the magnetic optics, proceeded in collaboration with numerous laboratories 
and universities (e.g., the Oak Ridge, Lawrence Berkeley, and Argonne national 
laboratories, the Northeastern University and the University of Texas) as well as 
with the industrial contractors associated to the project. In particular, a state-of- 
the-art RFQ was built at Los Alamos and delivered to serve as a proton source in 
the Superconducting Super Collider (SSC) expected to be built in Texas (103]{704]. 
Similarly, development and testing of the magnetic optics components were done 
at the Argonne and Lawrence Berkeley laboratories. 


However, although all its components had been built, and many of them suc- 
cessfully tested, the GTA was never fully assembled, and its construction aban- 
doned at the end of 1993 [104]. The floor space was later used for assembling and 
testing the Low Energy Demonstration Accelerator (LEDA) of the Accelerator 
Production of Tritium (APT) project [105], various equipments and spare parts 
were reused in other projects, and finally the GTA accelerator itself was donated 
in 2004 to the University of Indiana for its Low Energy Neutron Source (LENS) 
facility, to be used as a training ground for scientists who will later work at the 
$1.5 billion neutron source to be completed in 2006 at Oak Ridge (106)... 


9.3 Emerging neutral beam technologies 


As recalled in the previous section, the years immediately following the collapse of 
the Soviet Union coincided with the abandonment of the most visible components 
of the U.S. national neutral particle beam research and development program, 
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the ISE experiment and the GTA accelerator. Nevertheless, this did not mean 
that interest in neutral particle beam was lost, neither that development stopped 
there. In particular, the numerous told and untold reasons for terminating an 
attempt to demonstrate the viability of a neutral particle beam system in outer- 
space were much more of a technical and political nature than of a fundamental 
one. With a successful proof of principle experiment such as BEAR, there was no 
really compelling scientific reason to make another space-based experiment which 
would not have added very much to the understanding of the underlying physical 
issues. In fact, problems such as developing an appropriate kW to MW class 
energy source for the whole system, a suitable cooling system for the accelerator, 
or a highly efficient beam neutralizer, would not have much benefited from a 
crash-program to put GTA into outer-space. 


The points to be stressed are therefore: 


1. That the accelerators and associated technologies for neutral beam weapon 
systems are so closely related to those of accelerators and detectors used for 
fundamental nuclear and elementary particle research that their development 
does not need to be done in military laboratories. And: 


2. That the specific characteristics of neutral particle beams for military appli- 
cations such as ballistic missile defense should be discussed in relation to the 
full range of technologies available for accomplishing similar decoy/weapon 
discrimination and/or target destruction objectives [74]. 


Indeed, the most important and almost unique advantages of high-energy par- 
ticles are their ability to penetrate deeply into any target and to interact strongly 
(that is by inducing nuclear reactions) with any substance. This implies that the 
potential damage can be considerable, and that even for very low-intensity beams 
the secondary particles emitted in the nuclear interactions with the materials in the 
targets provide a signal giving a lot of information on their composition. This is 
illustrated by the current method which uses background cosmic-ray particles as 
a natural beam for the remote analysis of the elemental composition of artificial 
satellites [107]. While this method requires times on the order of months or years 
to integrate sufficient data to achieve useful data, directing a modest beam of ac- 
celerated particles at a given spacecraft would achieve the same result in a fraction 
of a second.* 

4This requires, of course, the use of very sensitive and clever techniques for detecting the 
secondary particles emitted by the target, and for discriminating them against the cosmic-ray 
background, something that can only be compared to the skill required to design the detectors 


and analyse the results of sophisticated high-sensitivity nuclear or elementary-particle physics 
experiments. 
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Considering that it is most probably easier and less expensive to destroy a 
ballistic missile or a spacecraft by means of some kind of a kinetic interceptor than 
by either a high-energy laser or particle beam, the most potent application of a 
space-based particle beam is quite certainly that of warhead/decoy discrimination. 
This is possibly why this application had been assigned in priority to space-based 
particle beam systems starting 1986 in the United States (see, e.g., references 
[98] and [99], where it was suggested that a discriminator/weapon decision was 
to be made by 1992). The remaining question to answer before designing and 
possibly deploying such discriminators is therefore whether there could be any 
strong competitor on the design horizon, of which two can easily be identified: 


1. Antimatter beams. While a neutral hydrogen beam containing 5-10 MeV 
protons would be sufficiently energetic to strongly react with the surface 
of a target, a beam with an energy of at least 50-100 MeV is necessary to 
penetrate inside, or to generate sufficiently many neutrons on the surface, 
in order to determine whether there could be a nuclear warhead within. If 
instead the beam would consists of low-energy antihydrogen atoms, the an- 
tiprotons would spontaneously annihilate on the surface, generating several 
high-energy pions for every antiproton hits, which would deeply penetrate 
into the target and strongly interact with the materials. An antihydrogen 
beam would therefore enormously simplify the design of the accelerator, 
which could operate at a much lower beam energy and current than for 
conventional hydrogen beam. Therefore, the technological burden would be 
transferred to the antimatter technology, which, however, is under intensive 
development since more that ten years [76]. 


2. Superlaser beams. Comparisons of the relative ability of neutral particle 
beam and laser systems for discriminating between reentry vehicle and de- 
coys show that particle beams can typically discriminate about hundred times 
as many objects as can lasers, and do so with significantly greater certainty 
[108]. This applies, however, to lasers with relatively low peak power. The 
recent invention of “chirped pulse amplification,’ which provided a factor of 
one million increase in the instantaneous power of lasers, enabled tabletop 
lasers to produce nuclear reactions directly [109]. Therefore, superlaser 
beams combine the ease in steering and focusing of optical laser beams with 
the capacity of particle beams to generate high-energy secondary particles 
in distant targets, which gives them the ability to “x-ray” remote objects and 
discriminate whether they are warheads or decoys [74. 


One can therefore conclude that neutral particle beams provide a credible 
option for discrimination, but that this option should be constantly compared to 
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emerging alternatives which arise as technology advances. In particular, it may 
happen that the laser system required to neutralize a negative hydrogen beam 
could in fact be of a complexity comparable to that of a superlaser able to do 
the same discrimination task on its own... But it would then be necessary to 
assess that further advances in antimatter technology would not require a complete 
reevaluation of the possible military uses of antimatter, both in offensive and in 
defensive weapons (see, [80] {96} (97) [76] and references therein). 


Conversely, if some major scientific or technical advance is made in accelerator 
technology, such as very-high-efficiency superconductive acceleration and radio- 
frequency generation, very-high-brightness laser-acceleration of neutral-particles, 
etc., neutral particle beams may find again their leading position as a potential 
rapid hard-kill system, simply because of the intrinsic strongly-interacting nature 
of high-energy particles.’ 


>The abstract of a typical comparative study (possibly biased in favor of particle beams) 
summarized this fact as follows: “This report explores the role of directed energy weapons 
(DEWSs) in theater defenses. For ranges shorter than 200-300 km they are much cheaper than 
space-based interceptors (SBIs); they are competitive with ground-based interceptors (GBIs). For 
inter-theater ranges of 1000km, lasers are competitive with the SBIs, but NPBs are significantly 
cheaper than either. For nominal laser and space-based interceptor (SBI) costs, lasers are strongly 
preferred for ranges under 300-500 km. For ranges 700 km, SBIs have a slight advantage. Neutral 
particle beams (NPBs) appear dominant for ranges over 400-1000 km” [110]. 
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Chapter 10 


Charged particle beam propagation 
experiments 


10.1 ATA at Lawrence Livermore National Labora- 
tory 


The purpose of the construction and operation of the Advanced Test Accelerator 
facility (the ATA) at the LLNL, as well as its basic characteristics, have been 
described in a number of informal, e.g., DIN (23) 112) OIZ], and more technical, 
e.g., (114) {115} [116], papers and reports. Summarizing from these publications: 


The main uncertainty in the concept of charge-particle beam weapons 
is whether it is feasible to propagate an intense self-focused electron 
beam through the atmosphere. That is, an electron beam held in tight 
focus by its own magnetic field. To conduct a comprehensive program 
of electron-beam propagation experiments, LLNL has constructed 
between 1978 and 1982 a 50 MeV, 10 kA, linear accelerator (the ATA) 
at its high-explosive test location, Site 300, which is well equipped for 
managing experiments with unusual hazards. 


Together with its associated program of beam propagation physics, 
the ATA represents the largest single component of the Defense Ad- 
vanced Research Agency (DARPA) particle-beam technology pro- 
gram, whose aim is to establish the feasibility of particle beam 
weapons. The prime goal of the Department of Defense particle- 
beam technology program is to resolve what is and is not possible in 
beam propagation. Accordingly, the goal of the ATA is to develop an 


71 


experimental capability that can resolve critical questions about beam 
propagation physics in a timely and cost-effective fashion. 


A first generation of particle-beam weapons will emphasize short- 
range applications. Potentially first applications of particle-beam 
weapons may be for the defense of large ships against cruise mis- 
siles. Another early use may be for terminal ballistic missile defense 
of hardened sites such as missile silos or national command authority 
centers. The short range, which reduces the sensor burden on search 
and fire-control radars, demands the high lethality that particle-beams 
possess. Boring its way out to targets at a rate of a kilometer per 
millisecond, the beam can deposit megajoule of energy almost instan- 
taneously. Consistent with fire-control system considerations, beam 
weapons can have the capability of engaging tens of targets per second. 
With such characteristics, charged-particle beams are particularly well 
suited to counter small, very fast, highly maneuverable threats. 


Photographs of the completed 200-meter-long ATA facility, as well 
as drawings of the 80-meter-long underground experimental tank are 
shown in reference [112]. Also visible is the 4-meter-thick shielded 
door which can be moved aside for beam experiments in open air. For 
such open air experiments the beam is directed towards a staging area 
where it may interact with various targets after propagating in free air 
over distances which are only limited by the topography of the ATA 
site, located in a shallow valley at Site 300. 


The main characteristics of the ATA in relation to his primary purpose, the 
study of endoatmospheric beam propagation and interaction with military targets, 
are as follows: 


Beam energy : 50 MeV 


An energy of 50 MeV means that that electron beam is fully relativistic but 
still substantially below the 1’000 to 10’000 MeV required for propagating a 
distance of one Nordsieck length in open air. This means that the experiments 
at ATA will be done under rather difficult conditions, which has the advantage 
that the results will easily extrapolate to higher energies. In this respect, it 
is important to recall the the main difficulties with particle acceleration is in 
the low energy section: as soon as an electron beam has reached an energy 
of 10-100 MeV it is easy to inject it into a betatron and further accelerate it 
to higher energies. An example of such a betatron, designed to accelerate 
the 10-kA ATA beam from 50 to 250 MeV, is given in reference [117]. 
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Beam current: 10 kA 


A current of 10 kA is on the order of the critical current given by equation 
(x), which corresponds to the maximum current for self-pinched propagation 
through the atmosphere a sea-level pressure. 


Pulse length : 70 ns 


A pulse length of 70 ns is adequate for studying nose erosion and tail losses 
during propagation, as well as the convective nature of major instabilities 
such the hose instability. 


Burst rate : 1 kHz — Average rate : 5 Hz 


A maximum repetition rate of 1°000 Herz, and an average rate of 5 Herz, 
imply that every second five successive pulses separated by a time-delay 
between | and 200 milliseconds can be sent into the experimental tank or 
into the atmosphere. This enables to study hole boring and stability of 
subsequent pulses through the low-density plasma generate by preceding 
pulses in reduced-density or ground-level air. 


Lethality : 35 kJ/pulse 


The is only 3.5 percent of the 1 MJ energy which would be delivered by a 
100 ns, 10 kA, 1 GeV beam pulse, i.e., the energy equivalent of about 0.25 
kg of TNT, which is considered to be on the order necessary to destroy a 
typical target. Nevertheless, this is sufficient to test the lethality of the beam 
on numerous targets without completely destroying them and the sensors 
used in the measurements. 


The accelerating principle used in the ATA is magnetic induction, a technology 
traditionally used in circular accelerators such as the betatron, which was pioneered 
and developed for use in linear accelerators at the Lawrence Livermore National 
Laboratory by N.C. Christofilos [118]. 


The construction of the ATA started in 1980. The first tests of the electron 
gun (injector) began in November of 1982, and its full 10 kA, 2.5 MeV, beam was 
delivered to the ATA main accelerator in January 1983. 


The ATA started operating at its full design specifications of 10 kA and 50 
MeV in July 1984. However, the beam current was was found to be strictly limited 
by instabilities within the accelerator so that only very uniform beam pulses could 
be accelerated to full energy when injected into the 85-m-long main accelerator 
structure. 
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The most serious such instability, termed beam break-up (BBU), is a very 
rapid growth of any beam’s transverse displacement to a disruptive amplitude. A 
radical cure to this problem was found and tested in the early 1985 using another 
lower-energy induction accelerator, the Experimental Test Accelerator (ETA), 
[120]. The idea was to fill the accelerator with benzene gas at a pressure 
of 10 * Pa, corresponding to an altitude of 80 to 120 kilometer above the Earth, 
and to created a plasma channel by sending a low-energy laser pulse through it. 
This proved to be very successful, the channel providing an electrostatic guide for 
the beam all the way through the accelerator [121]. Moreover, this breakthrough 
provided a way to greatly simplify the construction, and to reduce the weight, of 
future linear induction accelerator. Indeed, quoting from an October 1985 review 
of the ATA progress: 


“Clearly, the laser guiding technology gives a tremendous improve- 
ment in accelerator performances as well as simplifies accelerator 
operator and future construction (i.e., no longer needed are transport 
solenoid or steering magnets)” p.3145]. 


The final confirmation of the full and reliable operation of ATA by the use of the 
laser guiding transport technique was given on September 1986, (123) {124]]. From 
then on the ATA could be used for what he had been built for: To investigate the 
feasibility of an intense charge-particle beam as an endoatmospheric point-defense 
weapon. It turned out that the answer to one basic question of principle, namely 
that the propagation of the beam over long distances in open air is possible, came 
very soon. 


Indeed, in his State of the Laboratory statement of 1987, the director of LLNL 
was able able to highlight: 


“A major accomplishment in the laboratory’s beam research program 
was the first demonstration of open-air propagation of an electron- 
beam. Using the Advanced Test Accelerator (ATA), researchers were 
able to tailor the electron beam to permit stable propagation in the 


open air” [125] p.3]. 


In the Beam Research section of the same annual report, the following details 
were given: 


“Recently we successfully transported a high-intensity, 50-MeV elec- 
tron beam from our Advanced Test Accelerator (ATA) into condition- 


ing and diagnostic cells and then into free air. This test of stable beam 
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propagation in free air is the first of its kind at this energy, level of 
beam current (5 to 10 kA), and pulse rate (1 Hz). We have studied 
carefully the effects of beam parameters upon stability in air and have 
begun measurements of ancillary phenomena that will be important 
in assessing the practicality of using high-intensity electron beams as 
tactical weapons” [126] p.54]. 


However, none of the experimental details concerning these propagation ex- 
periments appear to have been published. Neither will there be any more detailed 
publications on the ATA and its uses over the following years until today. In partic- 
ular, there is no information whether the ATA beam has been used as an input into 
another, most probably circular, accelerator to increase its beam energy from 50 
MeV to 500 MeV or more [117] {127 [128]. The only subsequent open publications 
related to the ATA are those concerned with its use as a driver for a free-electron 
laser, e.g., [129], or publications related to technological developments which may 
be used to upgrade the ATA or to build highly reliable and efficient components 
for a new generation of high-current high-energy accelerators, e.g., [130]. 


In conclusion, if take the above statements for granted, we have to assume 
that the ATA project has succeeded in meeting its stated goals. This means, in 
particular, that single pulses must have propagated in a stable manner in free air 
over a distance of at least one Nordsieck length, equation (4.25), i.e., about 20 
meters according to the calculations made at the end of Section 4.3. Since ATA is 
able to fire five closely spaced pulses in a single burst, operation in this burst-mode 
must have allowed to verify that the stable propagation of a tightly focused beam 
is possible up to a distance of about 100 m. According to the available drawings 
of the ATA facility, e.g., the cover page of the brochure [LL], this range is about 
the distance between the exit port of the ATA and the focal point of the out-doors 
staging area. 


10.2 RADLAC at Sandia National Laboratory 


The main advantages of linear induction accelerators such as the ATA used for 
beam research, or such as the FXR, DARHT, or AIRIX used for flash x-ray 
radiography, are their intrinsic simplicity and capability to produce high quality 
beam pulses under reliable conditions. However, while such accelerators are well 
suited for research applications, they are very heavy and bulky, and therefore not 
suitable for applications in which relatively compact and light-weight accelerators 
are required. 
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One alternative technology has been successfully developed by A.I. Pavlovskii 
in the Soviet Union, with possible applications to linear induction accelerator 
without iron [72], and pulsed air-cored betatrons [132]. This technology, 
which does not use ferrite- or iron-loaded cavities, was subsequently developed 
in the United States, and the first device based on this principle was built and 
successfully operated at the Sandia National Laboratory [133]. The current and 
energy achieved in this accelerator, called RADLAC-I, were of 25 kA and 9 MeV 
with an average accelerating gradient of 3 MV/m. This accelerating gradient is 
substantially larger than that of the ATA (about 0.5 MV/m) which implies that 
the RADLAC is a promising candidate for a compact high-power accelerator. 
Moreover, the RADLAC can operate at high repetition rates since it is not affected 
by the classical beam breakup (BBU) instability, so that laser guiding of the beam 
in a laser-generated channel within the accelerator is not necessary p.1185].! 


The RADLAC technology was substantially improved during the 1980s and the 
construction of RADLAC-II initiated. This accelerator consists of two successives 
accelerating modules called RIIM. Using such modules in succession a beam can 
be accelerated to higher and higher energies. In 1985 RIIM was capable of reliable 
operation at output levels of 40 kA and 9 MeV, [135]. 


With the RIIM operating, RADLAC-II could be assembled and beam propa- 
gation experiments using its 40 kA, 18 MeV beam were soon successful [136]. In 
early 1986 the beam was extracted, without significant losses, and propagated into 
a magnetic-field-free, air-filled experimental tank. At a pressure of 1 atm the beam 
propagated straight without oscillations, and the radius was measured to be about 
0.75 cm, somewhat smaller than the 0.9 cm beam radius within the accelerator. 
The RADLAC-II beam was then conditioned in a 16-m-long, ion-focused region, 
and was allowed to propagate outside the accelerator building in open air where it 
propagated in a stable manner for quite a distance. As a matter of fact, figure 9 of 
reference is an open-shutter photograph of the RADLAC-II beam propagat- 
ing outdoors at night. The Manzano mountains of New Mexico are visible in the 
background. 


Subsequent publications gave no further details on out-doors propagation ex- 
periments. Nevertheless, research and development did not stop there, as is 
indicated by the summary presented at the subsequent DARPA conference [137] : 


“RADLAC program activities are reviewed. The work is broadly 
categorized under lead pulse stability (LPS), channel tracking, and 
Recirculating Linear Accelerator (RLA) activities. In LPS activities, 


'For a general discussion of this technology, and its applications to both linear and circular 
accelerators in the United States and in Russia, see the review [134]. 
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stable, open-air propagation of the RADLAC-II beam was demon- 
strated over ranges longer than a Nordsieck length. These shots were 
coordinated with the activities of other experimenters measuring beam 
induced emissions, and demonstrated that RADLAC-II could be fired 
on a predetermined schedule to allow numerous, coordinated, and 
geographically widespread measurements to be made. Since those 
experiments, improvements in the RADLAC-II accelerator, ion focus 
regime (IFR) beam conditioning cells, and matching of the accelerator 
beam to those cells have produced a beam which should allow greater 
than 20 betatron wavelengths in a Nordsieck length and saturation of 
hose growth to be observed. Channel tracking activities have included 
continued hardware development on the RADLAC-II Module (RIIM) 
for pulse-to-pulse channel tracking, the design of a laser for conduc- 
tivity channel tracking, and demonstration of a crude beam director for 
a high current beam. Codes which allow channel tracking simulations 
to be done have also been developed. Pulsed power and beam trans- 
port experiments on the Recirculating Linac have led to hardware and 
techniques which will allow demonstration of beam recirculation of a 
high current beam this year and a recirculating linear accelerator next 
year. These transport schemes and pulsed power developments can be 
extended to higher energies and a conceptual RLA for Navy charged 
particle beam weapon (CPBW) applications has been developed.” 


As well as at the 1987 SDIO/DARPA Services Annual Propagation Review 
[138] : 


“The RADLAC program encompasses high power electron beam 
propagation experiments and accelerator development, both for ad- 
vanced propagation experiments and to develop compact accelerator 
options for future charged particle beam weapons (CPBW). Propa- 
gation experiments include conditioning cell and lead pulse stability 
(LPS) experiments on RADLAC-II, and channel-tracking experiments 
on IBEX. The RIIM accelerator was used for two-pulse accelerator ex- 
periments to explore two-pulse configurations for RADLAC-II. The 
ion focused regime (IFR) transported, recirculating linear accelera- 
tor (RLA) experiment is aimed at future CPBW compact accelerator 
development. This paper briefly outlines recent work in these areas.” 


Concerning the stable open-air propagation of the RADLAC-II beam, it is 
remarquable that it is was demonstrated in 1986 already, [136], that is possibly 
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one year earlier than the same experiment with the ATA beam, [125] p.3] and (126) 
p.54]. This illustrates both the maturity of the RADLAC technology, and the level 
of the inter-US competition between national laboratories. 


As for the significance of the RADLAC-II open-air experiments in comparison 
to the corresponding ATA experiments, it should be remembered that what matters 
most in first order is beam power, as is shown by the elementary solution (4.24) 
to the Nordsieck equation (4.22). For ATA the initial beam power is Po 
05TW, and for RADLAC-II P) 0:72 TW. Therefore, according to Nordsieck’s 
approximation, (4.25), RADLAC-II should have a single-pulse range on the order 
of about 30 m, instead of about 20 m for ATA. However, if the calculations are 
made by computer-integrating the complete beam envelope equation using 
Moliére’s theory of multiple scattering and detailed energy-loss models, it 
turns out that both RADLAC-II and ATA have about the same efold range of 
approximately 22 m. In other words, while Nordsieck’s equation provides a 
good first approximation (especially for beams of high-energy particles), computer 
simulations are indispensable in the relatively low-energy domain in which both 
ATA and RADLAC -II are operating. 


A final important conclusion deriving from the RADLAC-II stable propagation 
experiments is that the plasma generation current (5.8), which should tend to limit 
stable propagation to currents less then about 10 kA, does not appear to be so 
critical, since the RADLAC-II beam intensity is of 40 kA. 


10.3 LIA-10 and LIA-30 at Arzamas-16 


As it originally started in the Soviet Union, e.g., (131) [132] [72], it is obvious that 
the development, and most certainly the use for beam propagation experiments, 
of iron-free linear induction accelerators continued in Russia. Indeed, the LIA-10 
accelerator of the 1970-1980s was upgraded to yield 50 kA for a 25 MeV pulse 
of 20 ns duration in 1993, and the construction of a new accelerator, LIA-30, 
producing a beam of 100 kA at 40 MeV, initiated [139]. 


However, there are no published results on long-range beam propagation ex- 
periments using these facilities. Moreover, while there is substancial published 
work on Russian beam propagation experiments through low-pressure air or even 
laser formed channels, e.g., [140 [141], there appears to be no open publications 
available on experiments similar to those made with ATA and RADLAC.? 


?The Rand Corporation report mentions that: “A follow-on report will discuss Soviet 
research on the propagation of intense relativistic electron beams through higher-pressure air and 
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10.4 PHERMEX at Los Alamos National Labora- 
tory 


A third technology suitable for making compact high-power accelerators is the 
more conventional standing-wave radio-frequency linac technology used in most 
high-energy linear accelerators built for fundamental nuclear and elementary- 
particles research. Using this technology, it is theoretically possible to accelerate 
a 10 kA beam pulse from 10 to 1000 MeV in a 30-meter-long linear accelerator 
[142]. 


This technology was in fact used in the first high-power flash x-ray facility, 
the PHERMEX accelerator complete in 1963, built at the Los Alamos National 
Laboratory to study the implosion of nuclear weapon’s primaries [143]. As this 
facility was expected to be superseeded by more powerful induction linacs (such 
as the DARHT, which is now operating) it was suggested in the early 1980s that 
PHERMEX could possibly be upgraded and used to study the endoatmospheric 
propagation of electron beams [22]. 


While the RF-linac technology has a few disadvantages compared to the 
induction-linac technology (e.g., a relatively short pulse length), it has an in- 
trinsic high repetition rate and multi-pulse capability. Besides, PHERMEX is 
located “in a blast-proof building at a remote, controlled access site [where] a 
clear line of site of approximately 2000 meters exists” [22] p.2]. 


In fact, some preliminary propagation experiments were performed using the 
available beam [144], and PHERMEX was subsequently upgraded to operate in 
the, up to 3 kA intensity, 20-60 MeV energy range which was anticipated to be 
theoretically possible (22) [145]. 


Emittance measurements for typical 300-500 A, 26 MeV, 3.3 ns micropulses 
were published [146], as well as several papers indicating that there are plans to 
inject the PHERMEX beam into a circular accelerator [128]. Besides 
increasing the beam energy, such an accelerator would have the advantage of 
providing a means to accumulate and condition the beam pulses before sending 
them into the atmosphere with a suitably larger energy, current, and duration. 


There are, however, no subsequent publications on these developments. In 
particular, it is not known whether or not PHERMEX is now used as an injector to 
a higher-energy circular (or even possibly linear, see p.2]) accelerator, and 
there is no published information on related beam propagation experiments. 


gases (P > 10 ? Torr)” 
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10.5 Proton and muon beams 


So far we have only considered experiments in which high-intensity electron beams 
are generated, accelerated, and injected into the atmosphere. The main reason for 
this was that electron-beam machines with currents in the kiloampere range still 
appear to be the only available ones to study endoatmospheric beam propagation 
over distances on the order of one or more Nordsieck lengths. Nevertheless, in 
view of the potentially larger range of proton versus electron beams, see Figures 
4.2 and 4.3, and other interests such as ion-driven thermonuclear fusion, research 
on high-current proton beam sources suitable for injecting a multi-kiloampere 
pulse into an accelerator is under way since quite some times. This is illustrated, 
for example, by the 6.0-kA, 1.1-MV proton beam obtained at Cornell University 
Laboratory of Plasma Studies [149]. 


Similarly, considerable research is underway for producing high-intensity low- 
emittance muon beams for fundamental research [150]. If the techniques used for 
producing these beams could be extrapolated to yield sufficiently powerful pulses 
of muons, a single such pulse could in principle propagate (according to Figure 
4.2 and 4.3) to ranges of up to several km in open air. The main problem with 
muons, however, is that they are short-lived particles, witha life-time of only 2.2 
microseconds at rest. Nevertheless, if the muons are accelerated, say to 1000 MeV 
(which since the mass of the muonism = 105 MeV/c? corresponds to a Lorentz 
factor of 10) they would be able to propagate to a distance of c 6:6 km 
before decaying. Therefore, in theory, a series of high-energy muon beam pulses 
could be used either to strike a target at a distance of tens of kilometers, or else 
to bore a channel in the atmosphere to guide a more powerful electron or proton 
beam to a distant target [76]. 


While such a concept may look very futuristic at present, it should nevertheless 
be seriously assessed in view of several synergistic factors which relate muon beam 
technology to other advanced beam technologies of importance to the subject of 
this report. For example, there is strong similarity between the techniques used 
to produce muons (see, e.g., chapters 2 to 4 in reference [150]) and those used to 
produce antiprotons (see, e.g., [80]) : In both cases the simplest method consists 
of striking a production-target with a sufficiently high-energy proton beam to 
generate copious amounts of either muons or antiprotons during the collisions, to 
collect as many as possible of these particles with a magnetic lens, and to decrease 
their kinetic energy dispersion (i.e, to “cool” them) in order to form a sufficiently 
low-emittance beam, which can then be stored or directly used for some purpose.’ 


3Another link between muon and antimatter technology, which was already mentioned in 
Chapter 8, is that antimatter is the sole portable source of muons, in the sense that upon annihilation 
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Moreover, if instead of producing the muons with a high-energy proton beam, 
a high-intensity electron beam is used instead, one is led to a system involving 
a multi-megawatt muon production-target driven by a multi-kiloampere electron 
beam, which requires the kind of technologies and experience associated with 
accelerators such as the ATA [I51].4 


In practice, the kind of muon production and cooling systems envisaged or un- 
der construction at present, e.g., [150] [153], are not expected to be fully operational 
before 2006-2010. Also, their muon yield will be orders of magnitudes below 
what would be needed to envisage concentrating a bunch of them into a pulse 
suitable for propagation in self-pinched mode over a significant distance in open 
air. Nevertheless, this situation is quite similar to that of antiproton “factories,” 
which have increased (and still increase) their output by many orders of magnitude 
over the years. 


every antiproton yields on average between 2 and 3 muons, which could be fed to a collecting and 
cooling system without needing a high-energy beam to produce the muons in the first place : 

4The first author of this paper headed the ATA project for a number of years, see [23] {112}. 
It should be remarked that in an electron-beam based system, the same high-power electron 
accelerator could be used as a driver for either a free-electron-laser or a muon-factory, which 
would both provide a means to guide the more powerful electron-beam towards a distant target. 
For more examples of such synergies, see, e.g., [152]. 
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Chapter 11 


Conclusion 


In this report we have reviewed the theory and the proof-of-principle experiments 
which have been carried out to demonstrate that high-intensity high-energy particle 
beam propagation in full or reduced-density atmosphere was possible, and that the 
range and stability of the beam pulses were in agreement with theory. 


We have, however, not analysed the practical consequencies of these conclu- 
sions, which determine with certainty the range of physical parameters compatible 
with realistic high-power directed-energy particle beam systems; neither have we 
investigated the R&D related to the possible construction and deployment of such 
systems. These are the subjects of a companion report,! of which only a small 
part, concerned with the radiological effects of directed high-intensity high-energy 
electron beams, has so far been published [154]. 


Nevertheless, by going through the theory and the details of several important 
experiments, we have followed several complete “research cycles,” that is several 
paths going from some initial ideas to their verification by means of a suitable 
experimental program. In the present case the outcomes of these research cycles 
can be qualified as “scientifically successful,’ meaning that the technical feasibility 
of some concepts have been established. While this is important, it does not mean 
that the associated technologies should necessarily be developed and applied, 
which is precisely the motivation for writing this report.” 


! Particle beam weapons: A review and assessment of current R&D, Report ISRI-82-05, to be 
published. 

“And let us not forget that a great breakthrough in military technology, like the invention of 
the H-bomb, can quickly come back to haunt us.” Hans Bethe, Physics Today, October 1978, p.13. 
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